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INVESTIGATION OF A NONLINEAR CONTROL SYSTEM 1 

By I. Flugge-Lotz, C. F. Taylor, and H. E. Lindberg 


SUMMARY 

Nonlinear elements are sometimes added to linear control 
systems in order to improve the response of the system to an 
arbitrary input. This can be done in different ways, one of 
them being the variation of the coefficients of the differential 
equation describing the system before the nonlinear elements are 
added. This variation of the coefficients may be done in a 
continuous or in a discontinuous way. In the present paper 
a discontinuous variation of the coefficients is studied in detail 
and investigated for practical use. 

The nonlinear feedback is applied to a second-order system. 
F rom former analytical considerations the process of control is 
visualized as establishing an ensemble of linear second-order 
differential equations ( some with stable and some with unstable 
homogeneous solutions) and switching from one equation to 
another so as to maintain small instantaneous error for rela- 
tively arbitrary inputs. Physically, this control process is 
realized with a linear second-order control system to which have 
been added possible discrete combinations of proportional and 
derivative feedback. The particular combination of feedback 
employed at any instant is determined by a feedback switching 
circuit which is in turn operated by sensed binary information 
obtained from the output, output derivative, error, and error 
derivative ( namely , the signs of these variables). Techniques 
that are common to the digital computer field are used to imple- 
ment this switching circuit. 

Once physical realization is completed, simulation techniques 
are used to study and evaluate the performance of the nonlinear 
control system and to compare it with a linear system for a wide 
variety of inputs. A detailed quantitative study of the influence 
of relay delays and of a transport delay is also given. In 
addition, the effects of physical imperfections that are likely to 
be encountered in any application of the control theory are con- 
sidered (e. g. , velocity and acceleration limits). 

An analysis of the experimental results shows that this type 
of nonlinear control system performs better than a linear control 
system having a natural frequency 15 times greater. For this 
comparison, performance is evaluated in terms of the average 
value of the magnitude of the instantaneous error for bawl- 
limited inputs. Further, in contrast with the linear system, the 
nonlinear system performance is virtually independent of vari- 
ation in the damping factor of the system. 

A preliminary extension of this type of nonlinear control 
concept to higher order systems is presented. Experimental 


results are given for a third-order system. These results show 
that just as in the second-order case the nonlinear system per- 
formance is better than that of a comparable linear system. 

INTRODUCTION 

With the demand for more exacting performance, more 
emphasis has been placed on nonlinear aspects of control 
systems. The term “control systems’’ can be interpreted 
to include active networks and feedback amplifiers as well as 
servomechanisms. From the standpoint of analysis, un- 
intentional nonlinearities have to be taken into account to 
explain performance. From the standpoint of synthesis, 
intentional nonlinearities have been introduced to improve 
performance. However, up to this date only in special 
cases have advancements been obtained in the field of 
nonlinear control systems. 

The designs of nonlinear control systems have inherent 
advantages. One advantage is that the response of a non- 
linear system at a certain time can be made less dependent 
upon past response than can a linear system of comparable 
power-handling capability. 2 This means that the non- 
linear system can be made to follow more arbitrary classes 
of inputs with less dynamic error than the comparable linear 
system. Another advantage is that the mathematical 
difficulties encountered may actually be conducive to con- 
sideration of more realistic criteria of performance. In the 
nonlinear realm it is essentially as easy to invoke a criterion 
such as the minimization of instantaneous error for non- 
stationary random inputs as it is to invoke the largest 
possible flat amplitude response for sinusoidal inputs. 

In the present paper a control system of second order, 
which was first suggested by Flugge-Lotz and Wuncli on 
the basis of analytical studies (ref. 1), is investigated. The 
physical realization of this system and its performance are 
studied in great detail by Dr. I. Flugge-Lotz and Dr. C. F. 
Taylor. A detailed quantitative study of the influence of 
relay delays and of a transport delay was made by Dr. H. E. 
Lindberg and is given in appendix A. 

Ibis in\ estigation was conducted at Stanford University 
under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. The 
authors wish to thank Dr. A. M. Peterson of the Electrical 
Engineering Department of Stanford University for his 
continued interest and his most helpful advice on the 


1 Supersedes NACA TN 3826, “Investigation of a Nonlinear Control System” by I. Flugge-Lotz and C F Taylor 1957 

2 In linear theory, the impulse response or the autocorrelation function of the system gives an indication of how pas’t response is weighed. 
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electronic problems which were encountered during this 
investigation. They also wish to thank Dr. G. S. Bahrs for 
his useful suggestions for a special transistor switching 
circuit. 


SYMBOLS 


A 

peak-to-peak amplitude of input 

a ,b,c 

parameters defining a system 

a + ,b + ,c + 

constants used in appendix B 

a,b,c, 

parameters in differential equation for control 
servo 

li 

viscous damping of motor and reflected load 
referred to motor shaft 

I) 

E=y—x=—e 

linear damping factor 

e 

instantaneous error, x—y 

i 

function 

( Ji, do 

input and output of Fade circuit, respectively 

H(j> ) 

transfer function, 1/(04?+ 1) 

/ 

inertia of motor rotor, gears, and reflected load 

i=V - 1 

k u k 2 

gain constants 

k m 

constant of proportionality between output 
velocity and back electromotive force 

M 

constant depending on initial conditions 

V 

operator, indicating differentiation with re- 
spect to independent variable; that is, d/dt if 
real time variable is used and d/dr if;iion- 
dimensional time variable is used 

T 

repetition rate or period 

t 

time 

td 

switching delay 


maximum allowable switching delay 

V 

input voltage 

X 

input into system 

y 

output from system 

y 

approximation of output 

a 

symbol used to denote different constants 

Pm = -iP sgn (?/e)-2/3 Sgn {y'e') 

(Pmylf n) min 

smallest values of parameters giving good 
nonlinear system performance 

1 18,2/3, 17, 27 

positive constants 

y n = — 17 sgn 

(ye)— 27 sgn (ye') 

b 

positive constant 

€ 

small positive quantity 


coordinates introduced in appendix E 

P 

radius of curvature 

a 

real part of complex frequency variable o-+?co 

7 

nondimensional time variable normalized with 
respect to a>„, c o v t 

ft 

nondimensional frequency, w/a>„ 

CO 

frequency 

0)y 

sgn (f)=f/\f\ 

natural frequency of undamped linear system 

( ) 

time average 

( r 

d( )/dr 

( y f 

d 2 ( )/dr 2 

0 

equality sign in equations which describe 


operations (see eqs. (2) and (3)) 


Subscripts: 

as 

d 

e 

i 

im 

L 

lin 

m, 71 = 0,1,2,3 

max 

min 

nonlin 

o 


at an actual switch point 
ideal or desirable 
error 

initial or input 

image 

limit 

linear 

maximum or upper bound 
minimum or lower bound 
nonlinear 
optimum, output 


REVIEW OF LINEAR CONTROL THEORY 

It is desirable to obtain from linear control theory some 
useful concepts that can be generalized to the nonlinear 
case. These concepts are: 

(1) Operational notation 

(a) Transfer functions 

(b) Block-diagram representations 

(2) Control criteria 

(3) Control through parameters 

At the onset second-order systems are considered. How- 
ever, there is no difficulty in extending these concepts to 
higher order systems. 


OPERATIONAL NOTATION. TRANSFER FUNCTIONS, AND BLOCK DIAGRAMS 

Consider a physical process or situation in which the 
output is described in terms of the input as 


d 2 y 

dt r 


a ,'i + b ^'1 -\-cy — x(t) 


( 1 ) 


where a, b, and c are constants, y=y(t ) is the output, and 
x=x(t) is the input. 

Utilizing the operator p=d , ,dt ) equation (1) may be 
written 

( ap 2 + bp + c) y = x (2) 

(Eq. (2) reads u (ap 2 +bp+c ) operating on y equals opera- 
tionally x”) Formal solution of equation (2) for the ratio 
of output over input yields by definition the operational 
transfer function for the system. Thus 


V-JL (3) 

x ap 2J rbp J rc 


The operational block diagram for the system is obtained by 
placing inside a box the transfer function, equation (3). 
Coming into the box is the input ; going out of the box is the 
output (see sketch a). 


X 

1 

Input 

ap 2 +bp+c 


V 

Output 


Sketch a 


Here transfer functions and block diagrams are utilized 
merely as shorthand operational notations for differential 
equations. This is opposed to the Laplace transformation 
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viewpoint where transfer functions (and thus block dia- 
grams) have the properties of functions of the complex 
frequency variable, p=a-\-iu. The reason for stressing this 
interpretation is that shorthand (operational) notation has 
proven useful in the transition to nonlinear control whereas 
the Laplace transformation viewpoint (e. g., synthesis in 
the complex frequency plane in terms of poles and zeros) 
has not. 

CONTROL CRITERIA 

To gage the performance of an actual system an ideal or 
desirable system is usually established as a straight through 
connection (i. e., yd=x as denoted in sketch b, a block dia- 
gram of an ideal system) . 


From the block diagram the differential equation for the 
open-loop system may be written: 


( 1 \ c hL+(JL\ 
\KJcJ dt^\K x kJ 



Similarly, the closed-loop differential equation is 


/ I \dhj (B+K 2 k m 
\KJcJ dt^\ K x k m 


)f+(i )y=x(t) 


( 4 ) 


(5) 


In either case the differential equation is of the form 

a (6) 


X 


Input 


Vd 

Output 


Sketch b 


Comparison between the desired output y d and actual out- 
put y is accomplished by utilizing the instantaneous error: 

e=(y<i-y) = (x-y) 

A control criterion or criterion of performance is defined as 
the minimization of some property of the instantaneous 
error e for a given class of inputs. The minimizing process 
can be exact (i. e., resulting from a variational formulation 
of the problem) or approximate. 

CONTROL THROUGH PARAMETERS 

In linear systems the process of control is usually physi- 
cally obtained by applying feedback and/or compensation 
to the system that is to be controlled. A control criterion 
is realized (as closely as possible) by adjustment of these 
applied quantities. The concept of control through param- 
eters is an interpretation of this control process in terms 
of the differential equation describing the process. A simple 
example illustrates this concept. 

Consider the position control servo shown in figure 1. 
The uncontrolled (open-loop) system consists of an amplifier, 
armature-controlled motor, gear train, and load. Closed- 
loop operation is obtained by utilizing proportional and 
derivative feedback. The gain constants K x and K 2 are 
adjustable. Armature inductance has been neglected. 



Figure 1. -Block diagram of simple positional servo. J, inertia of 
motor rotor, gears, and reflected load; B, viscous damping of motor 
and reflected load referred to motor shaft; k m , constant of propor- 
tionality between output velocity and back electromotive force (it 
includes armature resistance and gear ratio from motor shaft to 
load). 


where a set of three parameters a , b, and c completely 
characterizes the system. It is possible, then, to view the 
process of control in terms of these parameters. One starts 
with a parameter set (a,6,c) defining the uncontrolled system. 
A control criterion yields an optimum parameter set 
(a,b,c) 0 . Control (feedback and/or compensation) is intro- 
duced ideally making it possible to adjust (a,b,c) to ( a,b,c ) 0 . 
In the above example the gain constants K x and K 2 afford 
this adjustment. 

This adjustment of the coefficients may be done in a con- 
tinuous or a discontinuous way. In reference 2 Schmid and 
Triplett have described an interesting and efficient way to 
vary the coefficients of a basically linear system continuously. 

NONLINEAR CONTROL 

TRANSITION TO NONLINEAR CONTROL 

In the preceding section it has been mentioned that the 
process of linear control of second-order systems may be 
visualized as the adjustment of the parameter set ( a,b,c ) to 
the set ( a,b,c ) 0 . The term optimum was used in the sense 
that some criterion of performance was approached as closely 
as possible. 

It still seems logical in the transition to nonlinear control to 
hypothesize control through parameters. The transition is 
obtained by allowing the parameters to become functions of 
the output y(t) and the input x(t) ; that is, 

a^>a(x,y) 

b^>b(x,y) 

c-^c(x,y) 

The mathematical description of the system is now 

a ( x ,y)^+l>(x,y) ) jt+c(x,y)y=x(t) (7) 

Without knowing the specific nature of the functions a(x,y), 
b(x,y), and c{x,y) it may be seen that equation (7) is a non- 
linear, inhomogeneous, and/or nonautonomous differential 
equation. Mathematically, little in general can be said 
about the solution of equation (7) given the function set 
a(x,y), b(x,y), and c(x,y). It seems, then, even more hopeless 
to attempt a synthesis problem which involves both finding 
the function set [a(x,y), b(x,y ), c(x,y )\ 0 for a specified control 
criterion and then physically realizing the system described 
mathematically. 
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NONLINEAR CONTROL THEORY 

One analytical attack on the nonlinear control-system- 
synthesis problem has been made by Fltigge-Lotz and Wunch 
(refs. 1 and 3 to 5). They suggested varying the coefficients 
a , b, and c , not continuously, but discontinuously. That 
means that for t\<Ct<Ct 2 there is one set of coefficients, for 
there is another set of coefficients, and so on. The 
different sets of coefficients are chosen in advance, but the 
times ti for change from one set to another are determined 
by the value and the decrease or increase of the deviation 
( x-y ). In other words, the system is linear in any interval 
1, but is nonlinear in the whole. The transition at 
any switching time t t occurs with continuous values of y(t) 
and dy/dt, but discontinuous values of (dhyjdtf). 



Phase-plane techniques were used for studying appropri- 
ate sets of coefficients and the appropriate dependence of the 
switching times on the deviations. 3 The authors of references 
1 and 3 to 5 succeeded in finding a switching rule which 
assures good performance in a multitude of cases. Their 
control system is mathematically described in the following 
way: 

^+2Z?(l+0 m ) ~+(l+7 n)y=x(r) (8) 


e instantaneous error, e=(ya — y) = (x—y) 
( Y=d( )/dT 


Actually, equation (8) is a normalized form of the control 
equation derived by Fltigge-Lotz and Wunch. However, 
the notation has been somewhat changed. See appendix B 
for a comparison of notations and the normalization involved. 

The subscript convention is 


10 + 20 


73 — 17+27 


02 — 10 20 
01= 10 + 20= 02 


72 = 17 — 27 
7i=— 17+27 =— 72 


(9) 


00= — 10 — 20= — 03 70=— l7 — 27=— 73 J 

PROPERTIES OF EQUATION (8) 

The following properties of equation (8) are noteworthy: 

(1) Equation (8) is a piecewise linear but overall nonlinear 
differential equation. 

(2) The parameters and y n are stepwise switching func- 
tions of time (their implicit variable). This property is 
illustrated in figure 2. 


(3) The time of switching and the particular combination 
of the parameters 0 m and y n employed at any instant are 
explicit functions of the output y and the input x. Spe- 
cifically, they are determined by quantized information 
derived from the output, output derivative, error, and error 
derivative, namely, the sign of the products sgn ( y'e ), sgn 
(y'e'), sgn (ye), and sgn (ye'). 

(4) There are basically 16 m,n subscript combinations and 
thus 16 & m ,y n parameter combinations. However, a detailed 
study of property (3) shows that only 8 are allowed. The 
allowed combinations may be 


or 


m=n= 0, 1, 2, 3 
m^n 


where 

x(t) input 

y(r) output 

1) linear damping factor (when 0 w =y n =O) 

r nondimension al time variable normalized with re- 
spect to co„; that is, r=co„£ 

natural frequency of undamped linear system, 
#=0m=7w=O 

0m= —10 sgn (y'e )— 2 0 sgn (y'e') 
m=0 } 1, 2, 3 

7 «= — 17 sgn (ye)— 2 y sgn (ye') 
n= 0, 1, 2, 3, 

i 0,20 ,i 7,27 positive constants 


sgn 


(f)=J-=f+ liov f> 0 

U ’ l/l \ —1 for/<0 


3 For details the reader is referred to references 1 and 3 to 5. Reference 1 contains the ideas 
but is so condensed that the inquisitive reader will find it useful to read references 3 to 5, of 
which reference 5 is probably the most accessible. Figure 3, p. 12, and figure 30, p. 70, of 
reference 5 will help in getting acquainted with the phase-plane trajectory of an output. Some 
of the original studies are described again later in the present paper when the performance 
of the system is discussed. 


m+7i=3 

The reason for the “exclusion principle ,, on coefficient 
combinations stems from a desire to obtain mirrored-image 
outputs for mirrored-image inputs; that is, 

y<m(r) = — y{r) 

is desired when 

x im (T) = — x(t) 

(5) The control criterion that was employed in obtaining 
the functional dependence of 0 W and y n was 

\y*—y\=\*-y\=V\<* 

where e is a small positive quantity. This criterion of main- 
tenance of small instantaneous error between desired output 
and actual output enabled reduction of equation (8) to the 
approximate autonomous differential equation (see refs, 
and 3 to 5) 

y"+2D(l+p m )y'+y n y=(x— y) 

~0 ( 10 ) 

Thus phase-plane techniques could be employed to find the 
functional dependence of /?,„ and y n . 
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(6) Once the linear damping factor D is fixed the process 
of control is obtained by switching parameters p m and y n . 
Equation (8) consists of an ensemble of eight (see property 
(4)) linear differential equations with constant coefficients. 
The process of control may be visualized as the switching 
from one member of the ensemble to another. This switch- 
ing is determined by quantized information derived from the 
input and output (see property (3)). From another point 
of view (consistent with the approximation described in 
property (5)), the output y is to be forced to satisfy two 
conditions simultaneously, that is, both sides of equation 
(10). This is approximately possible by switching to 
various p m and y n parameter combinations and can be 
visualized as the process of switching to various phase 
trajectories of equation (10) in the phase plane of y' against y . 

DISCUSSION OF NONLINEAR CONTROL 

In the section entitled “Transition to Nonlinear Control” 
a logical transition to nonlinear control systems utilizing the 
concept of control through parameters is suggested. How- 
ever, mathematical difficulty hampers the development of 
this approach. In the next two sections a particular non- 
linear control theory is presented. This theory constitutes 
the first step in the synthesis of a nonlinear control system 
which obtains control through parameters. Since the func- 
tional dependence of the parameters has been established, 
the problem is reduced to finding a set of five (constant) 
parameters (D^P^P, i7, 2 7) 0 - It should be appreciated, how- 
ever, that even the optimization of this five-parameter set 
cannot in general be accomplished analytically because of 
the overall nonlinear nature of the problem. 

Aside from questions on the analytical optimization of 
parameters in equation (8), there are equally important 
practical questions such as: 

(1) Can a useful control system that is described by the 

nonlinear differential equation (eq. (8)) be 
realized? 

(2) If the system is realizable, what is its physical 

nature? 

(3) If the system is realizable, how does it compare in 

performance and complexity with a “good” 
second-order linear control system? 

There are then mathematical difficulties on the one hand 
and physical difficulties on the other. The mathematical 
difficulties could be handled by numerical methods of inte- 
gration of the differential equation (e. g., utilizing a digital 
computer). However, this would give little insight into the 
nature of a system that is controlled through discontinuous 
variation of the parameters p m and y n . It has been found 
advantageous to investigate the physical questions first and 
then to utilize simulation techniques (analog computer) to 
investigate the analytical properties of this type of control. 

PHYSICAL REALIZATION 
PHYSICAL MODEL 

It is desirable to study the nature of a physical control 
system that is described by an ensemble of eight linear 
differential equations with control being accomplished by 


switching from one member of the ensemble to another. 
To do this, equation (8) is rearranged as shown below: 

^+2Z»^+2/=x(r)-(2Z)/3 m ^+Tn2/) (11) 

or in operational notation 

(p 2 +2Dp+l)y=x—(2Dp m p+y n )y (12) 

Forgetting for the time being that the parameters p m and 
y n are actually functions, one interpretation of this opera- 
tional equation and thus of equation (8) is shown in figure 3. 



Figure 3. — Block diagram of equation (8) assuming that (3 m and y n 
are constant (denoted by encircling dotted line). 

Figure 3 can be modified to take into account the fact 
that p m and y n are stepwise switching functions of time, 
their implicit variable (i. e., p m and y n can each take on four 
discrete values). This is shown symbolically in figure 4. 
The explicit functional dependence of the parameters p m and 
y n has not yet been given and is thus indicated as a switch- 
ing logic of undefined character. 



Figure 4. — Block diagram of equation (8) taking into account 
stepwise switching nature of p m and y n . 
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Utilizing the block diagram of figure 4, the physical inter- 
pretation of the nonlinear control system described by equa- 
tion (8) is quite straightforward. This system consists of: 

(1) A linear feedforward portion. This portion could be 
a linear control system in itself (e. g., the simple position 
control servo of figure 1 and equation (5)). 

(2) A feedback switching circuit comprised of: 

(a) Four discrete values of proportional feedback y n 
(two positive and two negative as shown in fig. 2) 

(b) Four discrete values of derivative feedback 2 Dfi m 
(two positive and two negative as shown in fig. 2) 

(c) A switching logic which at any instant determines 
the particular combination of derivative and propor- 
tional feedback 2Dp m ,y n employed 


SWITCHING LOGIC 


In this section it is shown that digital-computer techniques 
can be utilized to establish the switching logic for the feed- 
back switching circuit mentioned in the previous section 
Recall that the parameters and y n have been defined 
as functions; that is, 

Pm= — 10 sgn (y'e) — 2 P sgn {y'e') " 


m= 0,1, 2, 3 

y n = — l y sgn {ye)— 2 y sgn {ye') 


y 


( 13 ) 


n= 0 ,l, 2,3 J 

where the subscript convention has been given by equations 
(9). Equations (13) determine the switching logic. Thus, 
for example, 0 3 is chosen when {y'e)<^ 0 and (y'e'X 0 and y 0 
is chosen when (ye)^> 0 and (ye')> 0, so that the combination 
0 3 Yo I s chosen when (i/d< 0, (|/V)<0, (ye)>0, and (ye')> 0. 
At this point it appears necessary to form the products 
ye , ye' , y'e , and y'e' and then to find the sign of these prod- 
ucts in order to establish the switching logic. Physically, 
however, the process of multiplication is to be avoided if 
possible. That there is a possibility of avoiding multiplica- 
tion may be gleaned by realizing that 

sgn (ab )= sgn (a) sgn ( b ) 


since 

ab __ a b 

wrwm 

Thus equations (13) may be rewritten as 

Pm =— sgn (y')[ 10 sgn {e)+ 2 psgn (e')]' 

m= 0 , 1 , 2, 3 

7n=— sgn (y )[ it sgn {e)+ 2 y sgn )e')] 

n= 0,1, 2, 3 


(14) 


7 o is chosen when 

y> 0, e>0, e'>0 

or when 

y< 0, *< 0 , ^'<0 

so that the combination (/3 3 ,y 0 ) is chosen when 
y> 0 , ?/< 0 , e>0 , e'>0 

or when 

y< 0 , y'> 0 , e<0 , e'<0 

From this example it can be seen that it is not necessary to 
find the signs of products but rather that it is sufficient to 
find separately the signs of y, y' , e, and e' . 

Since the sign of a variable is quantized binary information 
of the variable, it is convenient to utilize digital-computer 
techniques to further the switching logic. This may be done 
as follows: 

Let the convention be adopted that y^> 0 be represented 
by 0 (binary zero), y<^ 0 be represented by 1 (binary one), 
and similarly for y', e, and e'. If the ordered sequence is now 
established as 

(y, v’ , e, e') 

four-digit binary logic may be employed to encode equation 
(14). In particular, a binary coded decimal may be used 
(see table 1). 

It was mentioned in property (4) of the section “Properties 
of Equation (8)” that not all of the 16 possible P m ,y n 
parameter combinations were allowed under their definition. 
This was termed an exclusion principle on the allowed coeffi- 
cient combinations. One of the advantages of the suggested 
binary coding scheme of table 1 is that this exclusion principle 
is built into the code. To understand this, consider the 
example of the combination p 3 y 0 given previously. In the 
code language fi 3 y 0 is chosen when the binary number 0100 
or 1011 occurs. What is implied by this example is that a 
binary number and its complement must be identical (i. e., 
0000 = 1111, 0111 = 1000) as far as the switching logic is 
concerned. Thus out of the 16 possible four-digit binary 
numbers only the first 8 are unique. That is, in counting 
from 0 to 7 in a binary coded decimal, if complements are 
included then so are the other 8 possibilities, 8 to 15 (see 
table 1). 

The allowed p m ,y n parameter combinations along with the 
encoded logic of table 1 are summarized in matrix form in 
table 2. Examples are given to illustrate the meaning of 
the table. In general the allowed subscript combinations 
are 

m = n = 0, 1 , 2, 3 
or 

m n 


m +7i = 3 

REALIZATION COMPLETED 


Again the subscript convention is defined by equations (9). 

Thus, utilizing equations (13), for example, is chosen 
when 

V’> 0, e<0, e'<0 

or when 


t/'<0, e>0 , c'>0 


Now that equation (14) has been successfully interpreted 
(encoded) in binary-logic form (table 2), t lie realization of a 
feedback switching circuit utilizing this encoded logic is a 
typical digital-computer switching-circuit problem. As is 
characteristic of any synthesis process there will, in general, 
be many ways to design this feedback switching circuit. 
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Figure 5.— Block diagram showing complete physical interpretation 
of equation (8). CD, zero-coincidence detectors. 


The block diagram of figure 5 shows one design that com- 
pletes the physical interpretation of equation (8) along the 
lines started in figures 3 and 4. 

In figure 5 the feedback switching circuit consists of: 

(1) The four discrete values of both derivative and propor- 
tional feedback 2D0 m ,y n 

(2) A relay switching circuit that connects the proper 
feedback combination 

(3) Zero-coincidence detectors CD that drive the banks of 
relays to one position or the other depending upon the signs 
of the sensed variables 

It should be noted that, depending upon the application, 
other forms of sign-sensing devices and other switching de- 
vices such as diodes, transistors, electronic switches, and/or 
magnetic amplifiers could be employed to obtain other reali- 
zations of equation (8). In any case the following properties 
are basic to any realization: 

(1) The signs of the four variables y, y ' , e , and e' are 
sensed. This may be thought of as the process of “reading 
in” the four-digit binary logic of table 1. 

(2) On the basis of the 2 4 possible binary decisions the re- 
quired feedback combinations 2D0 mj y n as defined in table 2 
are connected around a linear second-order member. 

It is important to stress that the only types of nonlinear 
operations required in the realization of this nonlinear con- 
trol system are switching-type operations. In addition, all 
the switching is to be performed in feedback paths, which 
means that the switching can be done at low electronic 
power levels. These practical features are definite design 
advantages. Thus, in summary, it can be said that this 
type of nonlinear control system is not only physically 
realizable but also practical from an instrumentation stand- 
point. 

EXPERIMENTAL VERIFICATION OF NONLINEAR CONTROL 

THEORY 

DISCUSSION OF SIMULATION TECHNIQUES 

Simulation techniques were chosen as an experimental 
mode of investigation of performance of the nonlinear con- 
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trol system realized from equation (8). These techniques 
offer the following advantages: 

(1) Proximity to the actual control system. This means 
that the same practical features with regard to instrumenta- 
tion (see the section entitled “Realization Completed”) are 
exploited to fullest advantage. Thus, just as in the actual 
system, the only nonlinear device required for the simulated 
model is a binary logic feedback switching circuit (see ap- 
pendix C for details). The linear portion of the system is 
simulated on an analog computer. Here the only operations 
required are summations, two integrations, and one dif- 
ferentiation. These are all operations which an analog 
computer does well. It can be said then that the accuracy 
to which the simulated model simulates equation (8) de- 
pends primarily upon the realized feedback switching cir- 
cuit. The most essential type of imperfection to be ex- 
pected in this switching circuit is time delay in switching. 
Exactly the same type of imperfection will be met in the 
physical control system. Thus there will be more nearly a 
one-to-one correspondence between the simulated model and 
the actual system than between either and equation (8). 

(2) Convenience in experimental investigation. In order 
to characterize the output y of the nonlinear system com- 
pletely, a set of five parameters 

iPy iPf 20, l7> 27) 

and the input x must be specified. In the performance 
evaluation of the system it is necessary to be able to vary 
these characterizing quantities conveniently. Simulation 
techniques allow this. 

PRESENTATION OF EXPERIMENTAL RESULTS 

Figures 6 to 17 present experimental results obtained from 
the simulation studies of equation (8). Briefly, the results 
are presented as follows: 

Figures 6 to 12 compare the responses (output y and error e) 
of the nonlinear system with that of a linear system for 
various classes of inputs x. (In comparing the linear and 
nonlinear responses it will be noted that there is not exact 
synchronism of events because, with the available experi- 
mental facilities, it was necessary to obtain the two responses 
separately.) The linear system utilized is that which con- 
stitutes the feedforward member of the nonlinear system 
corresponding to the case where 0 m =7n=O. The nonlinear 
system for figures 6 to 15 is 

03= 0o = 2 

02= -ft = 0.5 

73=— 70=2 

72=— 71 = 0.5 

Figure 6 compares the system responses to sinusoidal 
inputs and figure 7 shows the responses to triangular-wave 
inputs. A partially integrated square wave x(t) is defined 
as the output of a first-order linear system characterized by 
the transfer function 77(p) = l/(ap + l), when the input 
X\ (t) is a square wave. The responses of linear and non- 
linear systems to this type of input are given in figure 8. 
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(a) Nonlinear system; ft = co/w„=0.1 to 0.4 in 0.1 steps. 

(b) Linear system; n = o>/w v =Q.\ to 0.4 in 0.1 steps. 

(c) Nonlinear system ; fi = co/w v =0.5 to 0.8 in 0.1 steps. 

(d) Linear system; f2 = o>/to»,= 0.5 to 0.8 in 0.1 steps. 

(e) Nonlinear system; ft = a>/w„=1.0 to 2.0 in 0.2 steps. 

(f) Linear system; ft = co/co,,= 1 .0 to 2.0 in 0.2 steps. 

Figure 6. — Linear and nonlinear system responses for 20-volt peak-to-peak sinusoidal inputs with frequency varied. D=0.6; 2.5 small divisions 
on time scale=l normalized time unit; tick marks at bottom of figures indicate when frequency was varied. 





INVESTIGATION OF A NONLINEAR CONTROL SYSTEM 


9 




(a) Nonlinear system; S2 = 2 tt/ 7 7 =0.1 to 0.4 in 0.1 steps. 

(b) Linear system; SI = 2tc/T=Q.\ to 0.4 in 0.1 steps. 

(c) Nonlinear system; ft = 27r/T=0.5 to 0.8 in 0.1 steps. 

(d) Linear system; U = 2tt/T=0.5 to 0.8 in 0.1 steps. 

(e) Nonlinear system; Q~2w/T=1.0 to 2.0 in 0.2 steps. 

(f) Linear system; ft = 27r/T=1.0 to 2.0 in 0.2 steps. 

Figure 7. — Linear and nonlinear system responses to 20-volt peak-to-peak triangular- wave inputs with period T varied. D= 0.6 ; 2.5 small divisions 
on time scale=l normalized time unit; tick marks at bottom of figures indicate when periods were varied. 
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(a) Nonlinear system; a=4 and 2. Tick mark at bottom of figure indicates when a was varied. 

(b) Linear system; a = 4 and 2. Tick mark at bottom of figure indicates when a was varied. 

(c) Nonlinear system; a=l and 0.5. Tick mark at bottom of figure indicates when a was varied. 

(d) Linear system; a = l and 0.5. Tick mark at bottom of figure indicates when a was varied. 

(e) Nonlinear and linear systems; a=0. Tick mark at bottom of figure indicates demarcation between nonlinear and linear systems. 

Figure 8.— Linear and nonlinear system responses to 20-volt peak-to-peak partially integrated square-wave input. D= 0.6; 2.5 small divisions 

on time scale=l normalized time unit. 


The responses of the systems to small sinusoidal inputs are 
given in figure 9. Figure 10 gives the responses to clipped 
sinusoidal inputs. Figure 11 shows the response of the 
nonlinear system to sinusoidal inputs that have been dis- 
placed with a direct-current component. Figure 12 shows 
the responses of linear and nonlinear systems to a triangular- 
wave input whose periods and amplitudes are randomly 
modulated. 

Figures 13 to 15 deal with the effects of imperfections that 
are likely to be encountered in the actual control system. 
Figure 13(a) gives the results and data of an experimental 
investigation on the effects of switching delays due to 


threshold in sensing the sign of the error, sgn ( e ) for a tri- 
angular-wave input. The experimental results for a con- 
stant 9-volt input (see fig. 13(b)) are given below: 

Threshold, mv 14 2G 36 44 

Peak-to-peak error, mv 44 124 220 290 

Figure 14 shows the effects of placing progressively smaller 
limits on the acceleration of the nonlinear control system. 
For each value of y" limit considered, the output, output 
derivative, output acceleration, and instantaneous error 
are shown. The effects of a velocity limit on performance 
of the nonlinear and linear systems are compared in figure 15. 
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(a) Nonlinear system; S2 = <o/co„=0.1 to 0.4 in 0.1 steps. 

(b) Linear system; 12 = co/w„= 0.1 to 0.4 in 0.1 steps. 

(c) Nonlinear system; 12 = a>/a>„ = 0.5 to 0.8 in 0.1 steps. 

(d) Linear system; 12 = co/co„=0.5 to 0.8 in 0.1 steps. 

(e) Nonlinear system; 12 = a>/co„= 1.0 to 2.0 in 0.2 steps. 

(f) Linear system; 12 = co/co„= 1.0 to 2.0 in 0.2 steps. 

Figure 9. — Linear and nonlinear system responses to 4-volt peak-to-peak sinusoidal inputs with frequency 12 varied. D = 0.6; 2.5 small division 

on time scale=l normalized time unit; tick marks at bottom of figures indicate when frequency was varied. 
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(a) Nonlinear system; £2 = w/a>„ = 0. 1 to 0.4 in 0.1 steps. 

(b) Linear system; £2 = co/w„=0.1 to 0.4 in 0.1 steps, 

(c) Nonlinear system; £2 = co/a>„=0.5 to 0.8 in 0.1 steps. 

(d) Linear system; £2 = co/co„=0.5 to 0.8 in 0.1 steps, 

(e) Nonlinear system; £2 = co/co„= 1 .0 to 2.0 in 0.2 steps, 

(f) Linear system; £2 = w/a>„= 1.0 to 2.0 in 0.2 steps. 


Figure 10. Linear and nonlinear system responses to 20-volt peak-to-peak sinusoidal inputs with frequency £2 varied that are clipped symmetri- 
cally at ±6 volts. D = 0.6; 2.5 small divisions on time scale= 1 normalized time unit; tick marks at bottom of figures indicate when frequency 
was varied. 
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(a) ft = co/co„=0.1 to 0.4 in 0.1 steps. 

(b) fi=co/co„= 0.5 to 0.8 in 0.1 steps. 

Figure 11. — Nonlinear system responses of 20-volt peak-to-peak 
sinusoidal input with frequency Q varied that has been displaced by 
— 10 volts (direct-current component) for obtaining osculation. 
1) = 0.0; 2.5 small divisions on time scale=l normalized time unit; 
tick marks at bottom of figures indicate when frequency was varied. 


Figure 16 presents the responses of some special eases 
of the nonlinear system (special with respect to the choice 
of the p m and y n parameter values). In this figure the 
system responses to a triangular-wave input are given for 
four different parameter sets. The parameters pertinent 
to these results are listed in table 3 for easy reference. 

Figure 17 gives the response of the nonlinear system pos- 
sessing a low linear damping factor D=0.1. The responses 
are for a triangular-wave input whose frequency was varied 
in the same manner as that of figures 7(a) and 7(c). 

Detailed discussions of these results are given in the 
section entitled “Discussion of Results. ” 

DISCUSSION OF RESULTS 

PERFORMANCE EVALUATION OF SINUSOIDAL INPUTS 

To complete the synthesis of the nonlinear control system 
that has been derived from equation (8), it is necessary to 
choose the magnitudes of the parameters that characterize 
the system, that is, (D fl p, 2 / 8 , 17 , 27 ) or (D,p s ,p 2 ,y 3 ,y 2 ). It is 
expected that the performance of the system depends on 




Figure 12. — Linear and nonlinear system responses for triangular- 
wave input whose period and amplitude are randomly modulated. 
D = 0.6; 1 small division on time scale=l normalized time unit. 




(a) Nonlinear system response to triangular-wave input with period 
varied as in figure 7 (a). D= 0.6; tick marks at bottom of figure 
indicate when period was varied. Threshold in sensing sgn e was 
44 millivolts. 

(b) Constant input x(t) was 9 volts. 

Figure 13. — Effects of switching delays due to threshold in sensing 

sgn e. 
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(a) Nonlinear system; no y" limit. 

(b) Nonlinear system; ± 12-volt y" limit. 

(c) Nonlinear system; ± 9-volt y" limit. 

(d) Nonlinear system; ±7£-volt y" limit. 

(e) Nonlinear system; ± 6-volt y" limit. 

(f) Nonlinear system; ±4|-volt y n limit. 

Figure 14. — Nonlinear system responses to 20- volt peak-to-peak sinusoidal input with varying acceleration limits. Z) = 0.0; n=a>/a>„ = 0.5; 

2.5 small divisions on time scale=l normalized time unit. 


the choice of these parameters. For studying their in- 
fluence, simulation techniques proved to be very convenient. 
Experimental results were given in figures 6 to 12 where 
the response (output y and error e) of the nonlinear system 
was compared with the response of a linear S 3 7 stem for a 
variety of inputs. This gives the possibility of establishing 
the properties of the nonlinear system not only by itself but 
also with respect to a linear standard. The linear system 
employed for this purpose was that which constitutes the 
ieedforward member of the nonlinear system (i.e., p m =y n =0). 


These experiments allow parameter values for good per- 
formance of the nonlinear system to be found. Analytical 
and practical considerations that aid in the optimization are 
treated later in the section entitled “Choice of Parameter 
Values. ” 

The sinusoidal responses of the nonlinear and linear systems 
are compared in figure 6. Here, the frequency range con- 
sidered was 0.1 ^ £2=o>/co, ^ 2. The peak-to-peak input ampli- 
tude was 20 volts. These results show that the nonlinear 
system reproduced the sinusoidal inputs up to the frequency 
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(a) Nonlinear system, 
(b) Linear system. 


Figure 15. — Linear and nonlinear system responses to 20-volt peak-to-peak sinusoidal input with ±4£-volt velocity limit. Z> = 0.(i; H=2tt/T= 

0.5; 2.5 small divisions on time scale=l normalized time unit. 


12=1.4 with virtually no instantaneous error when com- 
pared with that of the linear system. For higher frequencies 
the error for the nonlinear system increased rapidly to the 
same order of magnitude as that of the linear system. 

Substantially the same comparative performance was 
displayed by the two systems when the input was a smaller 
(4-volt peak-to-peak) sinusoid as is shown in figure 9. 

Figure 11 gives the response of the nonlinear system to a 
20-volt peak-to-peak osculating sinusoid over the frequency 
range 0.1^12^0.8. This is a severe type of input for the 
nonlinear system since both x and x' simultaneously go to 
zero. This implies that y and y' are also small so that in 
equation (8) the discontinuous variations of the parameters 
p m ,y n cannot be so effective in determining the acceleration 
y", since 

y''=x-[2D{l+fi m )y' + (l+y u )y] (15) 

From the figures it is seen that the nonlinear system did have 
some difficulties near the osculating regions; further, the 
peak error increased as the frequency was increased (see 
appendix D). As might be expected from equation (15) 
this error for small values of y and y' can be reduced by 
increasing the magnitudes of the parameters (this will be 
discussed in more detail in the section entitled “Choice of 
Parameters”)- At any rate, by comparison, the error for 
the present system is always less than that of the linear 
system of figure 6. (In the linear realm, error is independent 
of a shift in the direct-current level of the input.) 

Taken collectively, the 20-volt, the 4-volt, and the oscu- 
lating 20-volt sinusoidal inputs tend to form a more realistic 
appraisal of the nonlinear system performance than does a 
single input amplitude. There is still much that can be 
learned from a detailed study of these three responses, but 
first it is desirable to obtain some sort of a quantitative com- 
parison between the nonlinear and the linear system per- 
formance. One such comparison can be made as follows: 

(1) Assume that the input to both systems x(t) is, and 
has been for a long time, a 4-volt peak-to-peak sinusoid of 
frequency 

£2=a)/a>„ = 0.1 

so that as far as the linear system is concerned this is a 
steady-state alternating-current input. 
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(2) Determine how much the band width or the natural 
frequency a>„ of the linear system must be increased in order 
that the time average of the magnitude of the instantaneous 
error 


= o r 2ir '° 

27rJo 


\e\dr 


for the linear system be reduced to that value given by the 
nonlinear system. 

Here it is easy to show that for the low-frequency steady- 
state alternating-current case 


f=T 4 D 


l^lu'n — ^ !^|/waz (15) 

7T 

For an input x=\x\ maz e inT the steady-state error is given by 

r ( 1 ~i+;a2Z)-a 2 )' 


or 

but 


e=x—y=\xl 

For small values of £2 this equation yields 

e~\x\ max [l-(l-iQ2D+ . . .)]e tQ * 
**\x\ max (itt2D)e iQT 

! ^ I max ~ 2T) 12 1 X I 


O^T 


therefore 


\*\itn=l r 

T Jo 


\e\max sin 12 t (/£2t j^lwaa; 

7T 


o|«l» 

7T 


(17) 


For the given system and input 



For the nonlinear system and the same input (approximate 
calculation from the much larger and clearer original of 
fig. 9 (a)) 


l^lnorciin ~ 0*91 

Thus, 

\e\nonlin~Y£ |^!//« 


Now from equation (17), since £2= co/co„, it is seen that in 
order to reduce \e\ un by this factor of 15 that co„ must be 
increased by the same factor. 
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(a) Complete system 1 and system with no derivative feedback 2. 
System 1: 0 3 = — 0o= 2; 0 2 = — 0i = O.5; 73 = — yo = 2 ; 72= — 71 = 0 . 0 . 
System 2: 0 3 = — 0 o =O; /S 2 = — 0 i = O; 73= — 7 o= 2 ; 72=— 7 i = 0.5. 

(b) System with no e' sensing 3 and system with no e sensing in deri- 
vative feedback loop and no e' sensing in proportional feedback 
loop 4. 

System 3; /S 3 = — /S 0 =2; /S 2 = — £i = 2; 73=— 70 = 2 ; 72=— 7 i = 2. 
System 4: 0 3 = — 0o = 2; /S 2 = — 0i = — 2; 73 — — 7 o— 2 ; 72 = — 71 = 2. 

(c) Complete system 5 and system with no derivative feedback 6. 
System 5: /S 3 = — /?o = 10; /S 2 = — 01=1; 73= — 70= 10; 72= — yi = 1. 
System 6: 0 3 = — 0o=O; 0 2 = — 0i = O; 73=— 70= 10; 72 = — 7l = 1. 

(d) System with no e' sensing 7 and system with no e sensing in deri- 
vative feedback loop and no e' sensing in proportional feedback 
loop 8. 

System 7 1 (3$= — 0o = lO; 0 2 = — 1 = 10; 73= — 70= 10; 72= — 71= 10. 
System 8: 0 3 = — 00=10; 0 2 = — /Sj = — 10; 73= — 70= 10; 72= — 71= 10. 

Figure 16. — Continued. 



In the example of the linear control system of figure 1 and 
equation (5) this increase in u> v by a factor of 15 would mean 
that the gain K x must be increased by the order of magnitude 
of (15) 2 since 

The gain K 2 must be increased by an order of magnitude of 
15 if the value of D shall not change. Such an increase in 
the loop gains is frequent^ not at all physically possible. 



Figure 16. — Responses of nonlinear system with four different para- 
meter sets to a 20-volt peak-to-peak triangular-wave input with 
fixed period T= IOtt. Z) = 0.6; 2.5 small divisions on time scale = 
1 normalized time unit; tick marks at bottom of figures indicate 
demarcation between systems. For systems’ numbering see table 3. 


(a) SI=2tt/T=0. 1 to 0.4 in 0.1 steps. 

(b) fl=27r/T=0.5 to 0.8 in 0.1 steps. 

Figure 17. — Nonlinear system responses to 20-volt peak-to-peak 
triangular- wave input with period 7 1 varied. Z> = 0.1; 0 3 = — 0o=12; 
0 2 = — /Si = 3; 73=— 7o=2; 72=— 7 i = 0.5; 2.5 small divisions on time 
scale =1 normalized time unit; tick marks at bottom of figures 
indicate when period was varied. 
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Up to this point little attention has been given to the 
detailed nature of the nonlinear system response. Closer 
inspection of, for example, figures 6 (a) and 6 (c) shows that 
the output y is a function that links serpentine fashion 
(oscillates) at a very high rate about the input x , but still 
the magnitude of the error is small. In fact, it is necessary 
to inspect the error at a scale 20 times larger than that of 
the output even to notice this phenomenon. Mathemati- 
cally this means that the functions x and y approach one 
another closely but that their derivatives differ appreciably. 
Physically, however, this is not at all undesirable as long as 
the magnitude of the error is small. (Actually for mechanical 
systems this property would be useful in preventing static 
friction.) This fine-grained oscillating character of y is the 
very essence of the nonlinear control theory. Every time 
the error or error derivative goes through zero the parameter 
set 0 m , 7n of equation (8) changes discretely as defined in 
equations (13) or by the binary logic of table 2. The discrete 
changes in the parameters cause discontinuities in the second 
derivative y " , which when integrated twice give y its serpen- 
tine character. To illustrate this point, the sketch of figure 
18 shows samples of the superposed input and output of the 
nonlinear system. The input in this case could be that of 
figure 6 (a) or 9 (a). 



Circle (a) in figure 18 is a typical cycle of the oscillating 
character of the output y. Commencing at r 0 , the error 
changes sign at r 0 , the binary number 0011 is “read into” 
the feedback switching circuit of figure 5, and using the 
notation of table 2 the parameter combination 0 373 is switched 
into the circuit. This causes an immediate reversal in the 
sign of the output acceleration y", so that at time r x the 
error derivative changes sign, the binary number 0010 is 
read in, and the parameter combination 0 272 is switched into 
the circuit; acceleration is still in the same direction but 
weaker. At time r 2 the error again changes sign, the binary 
number 0000 occurs, and the combination 0 O 7o is switched 
into the circuit; acceleration is in the opposite direction. 
At time r 3 the error derivative changes sign, the binary 


number 0001 occurs, and 0 l7l is switched in; this reduces 
the acceleration until at time r 4 the error again changes 
sign, 0011 occurs, 0 373 is again switched in; the cycle is 
complete. Although it was not mentioned at the time, 
figure 2 shows this sequence of 0 m , 7n combinations. The 
nature of the acceleration resulting from the switching can 
be seen in figure 14 (a). 

The comparable switching cycle for ?/>0 and ?/'< 0 as 
shown in circle (b) of figure 18 is: 

Logic Parameters 

0111 A)73 

0110 0 l72 

0100 037o 

0101 0 27l 

0111 0O 7 3 

Similarly, for ?/< 0 and y'<^ 0 as shown in circle (c) of figure 

18: 


Logic 

Parameters 

1100 

0373 

1101 

0272 

mi 

0o7o 

1110 

0i7i 

1100 

0373 


This is the mirror-image switching cycle for circle (a) (see 
property (4) in the section “Properties of Equation (8)” and 
also the discussion in the section “Switching Logic”). 

Finally, for y <^ 0 and y '^> 0 as in circle (d) of figure 18: 

Logic Parameters 

1000 0 o73 

1001 0 l72 

1011 0 37o 

1010 0 27l 

1000 0 o7 , 

This is the mirror-image switching cycle for circle (b). 

With this insight into the detailed behavior of the non- 
linear system more information can be obtained from the 
experimental sinusoidal responses of figures 6, 9, and 11 that 
have up until now been treated from a macroscopic rather 
than microscopic viewpoint. Along these lines, the following 
experimentally observed facts are noteworthy: 

(1) In comparing the errors for the 4-volt and the 20-volt 
peak-to-peak sinusoidal inputs in the frequency range of 
good reproduction 0.1^12^0.8, it is seen that: 

(a) The magnitudes of the errors are nearly the same (see 
figs. 6 (a) and 6 (c) and 9 (a) and 9 (c)). 

(b) The period of the error is generally smaller in the 
larger amplitude case. This is even more pronounced if the 
4-volt peak-to-peak case is compared with figure 11 for y 
in the region of —20 volts. 

(2) Good reproduction is characterized by many e and e ' 
switchings (see fig. 18) per cycle of the input. As input 
frequency is increased and the upper limit of small error is 
reached (fig. 6 (e) and 9 (e)) the e and e ' switchings become 
more infrequent until there are finally only two of each per 
cycle of the input. 
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The fact that the period of the error is smaller for larger 
inputs can be gleaned from equation (8) when it is rearranged 
as 

— [ 2/^(1 Pm)y f (# — y ) 

or roughly 

y" ~ — [2Z)(1 -\-p m )x' -\-y n x] (18) 

For a system with specified P„ n y n parameter values, the larger 
input and input derivatives will give stronger discontinuities 
in y" as the / 3 m ,yn combinations change. This implies that 
Pm,y n will change more often making the period of the error 
smaller. See figure 14 (a) for the justification of the approxi- 
mations in equation (18) since ?/, y ' , y " , and x are shown in 
this figure. (The input x was sketched in by hand.) 

In the frequency range of good reproduction the reason 
that the magnitude of the error is virtually independent of 
the input amplitude cannot be explained from equation (8) 
since this fact is intimately tied in with the imperfections 
in the feedback switching circuit (see section entitled 
“Effects of Switching Imperfections” for details). As 
measured from the experimental sinusoidal responses, 
\e\ max is of the order of 15 to 20 millivolts. It might be 
noted that earlier in this section the smaller 4-volt peak-to- 
peak input sinusoid was employed in the comparison when 
it was determined that co, should be increased by a factor 

of 15 in order to obtain the same \e\u n . This choice of the 
smaller input was decidedly in favor of the linear S 3 7 stem 
since the linear error increases linearly with input amplitude. 
Thus, if the 20-volt peak-to-peak input were used, an 
increase in co v by a factor of 5X15 = 75 would be required 

to obtain the same \e\u n . 

The fact that there is a correlation between good repro- 
duction (small error) and the existence of many error and 
error-derivative switchings is very important since it is the 
key to understanding the upper limits of good performance 
of the nonlinear system. One approximate way to investi- 
gate this matter is to establish a deficiency between the 
output acceleration required for good reproduction (many 
switchings) and the available output acceleration. Equation 
(18) gives an approximation to the output acceleration when 
the system is functioning well. As an example of this 
approach consider the input x to be a sinusoid of frequency 
SI (fig. 6 or 9). Since 

(•E | max = ^ 

it is to be expected that the nonlinear system will have the 
greatest difficulty in the vicinity of \x\ max where, from equation 


Thus, if in this region the output y is going to interweave 
the input x (as is characteristic of good reproduction), then 

\V | max 

This inequality then places an upper limit on good perform- 
ance of the system in response to sinusoids 

o<VW max (19) 


For the system used in obtaining the experimental results of 
figures 6 and 9, \y n \max =< 2) therefore, 

a<V2 

Inspection of both figures 6 (e) and 9 (e) tends to substantiate 
the above result. For example, in figure 9 (e) it is seen for 
11=1.2 (after the transient caused by turning on the input 
has been absorbed) that as the input goes through its maxi- 
mum, frequent e and e' switchings stop and do not occur 
again as the frequency is increased. 

PERFORMANCE EVALUATION OF INPUTS OTHER THAN SINUSOIDAL 

In the preceding section only the sinusoidal response of 
the nonlinear control system was discussed in studying 
performance. However, it is easy to see that the nonlinear 
system will cope with any input in the same manner as it 
does with sinusoids. Thus, as long as the magnitude of 
input acceleration does not continuously exceed the available 
magnitude of the output acceleration the switching process 
will commence and excellent reproduction will result. This 
type includes inputs with discontinuous derivatives and 
discontinuous inputs. The experimental sinusoidal re- 
sponses of figures 6, 9, and 11 themselves give some indica- 
tion of these facts. For example, in figure 6 (a) there was 
an initial discontinuity in the input and there were dis- 
continuities in the input derivative when the frequency 
was changed. The results in figures 7, 8, 10, and 12 prove 
further that the nonlinear system response is not dependent 
upon any specific type of input. Given in these figures are 
triangular-wave, partially integrated square-wave (including 
square-wave), clipped sinusoidal, and random 4 inputs, 
respectively. 

It should be noted that in the literature (refs. 1 and 3 to 5) 
a marginal-type input that would present a case of indecision 
to the switching circuit of the nonlinear system is discussed. 
This case has never been encountered experimentally even 
when the attempt was to produce this case. Thus, the 
marginal-type input is not considered practically important. 

USE OF PHASE-PLANE METHODS TO STUDY PERFORMANCE 

The phase-plane methods that were used in the original 
analytical development of equation (8) (refs. 1 and 3 to 5) 
can also be gainfully employed in studying the performance 
of the nonlinear control system derived from equation (8) 
once the and y n parameter values have been specified. 
(In this section it is important to distinguish between jS m 
and y n parameter values, i. e., /3 3 =2, /3 2 = 1, etc., and @ m y n 
parameter combinations, i. e., /3 3 y 0 , etc.) 


< In the strict statistical sense the probability that the inputs shown in figures 12 (a) and 
12 (b) are not samples of random stationary time series is admittedly high because of existing 
experimental facilities. These inputs were obtained by random manual modulation (both 
frequency and amplitude) of a triangular-wave input of peak-to-peak amplitude A and period 
T where 

0^A<!20 volts 
0.1 ^fi=27r/ 7^0.8 

Random manual modulation means that the operator varied by hand both the frequency 
and amplitude controls of the input generator as randomly as possible. In the present in- 
vestigation the comparative results of the linear and nonlinear responses to what appear to 
be band-limited random inputs are felt to be more important than the exact statistical prop- 
erties of the inputs. 
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If in equation (8) the error e=(x—y) is assumed small, 
then the output y can be approximated by portions of 
curves that satisfy the autonomous differential equation 

042Z?(1+|BJ^-+Y.y=o (20) 

where D , /3 m , and y n are defined in the section “Nonlinear 
Control Theory.” Further, if dyldr=y', these approxi- 
mating curves are defined in the y'y phase plane by 

(y*2-y'P=M(y\-y'f' (21) 

where M is a constant depending on initial conditions and 

Xl,2='~^(l ±V^ 2 (l+^m) 2 — 7 n 

Equation (21) comes from integration of the first-order 
differential equation 

dy' = — 219(1 +P n ) y '—y n y 
dy y\ 

= -2D(l+p m )-y n lj 7 (22) 

Since the particular f$ m , y n combination employed in each 
point of the phase plane depends upon y, y', e, and e', four 
approximating curves go through each point (see appendix 
E). The tangents to these curves (eq. (22)) indicate four 
directions which lie in an angular sector (see fig. 19). This 
angular sector is defined by the two extreme directions which 
apply if e and e' have the same sign. The two inner direction 
values apply if e and e f have opposite signs. By superposing 
the input x'x phase plane on this y'y plane it can be stated 
that at any point the tangent to the input phase curve must 



Figure 19.— Phase-plane angular sectors defined by equation (22) for 
the given parameter values. i/3=1.5; 2 0 = O.5; 17 = 0.45; 2 y = 0.05; 
2D = 0.25; w„=2. 


be included within the shaded angular sector (discussed 
above) if small error is to be obtained. Thus, equation (22) 
can be used to study performance. The width of the angular 
sector changes with the ratio y/y' for given and y n param- 

eter values. It is largest (180°) for y \— >0 and smallest for 
y = 0 if (1+j 3 m ) and y n have the same sign. However, the 
width for y '-> 0 is not significant since all phase curves 
(including the input) have infinite slope there. Thus, along 
the line y' = 0 the curvature of the input phase curve must be 
used to determine limits on good performance. The radius 
of curvature of the input must be smaller than that of the 
flattest approximating curve. 

In order to illustrate the use of these phase-plane methods 
in predicting limits of good performance consider the example 
of a sinusoidal input 2 = sin ftr. In this case the input is 
represented by an ellipse in the phase plane and 

<*(»') x 
dx “ x' 

In the superposed xx' and yy' phase planes of figure 20 are 
shown three sinusoidal-input phase curves (i. e., three 
different frequencies) and the families of phase trajectories 
of equation (22) for the 0 m , y n combinations where e and e' 
have the same sign. (In order to avoid extensive computa- 
tions the parameter values indicated in fig. 20 are those 
of an earlier investigation (see ref. 1).) Thus, the tangents to 
two intersecting phase curves define the angular sector at that 
point, as has been discussed. By tracing the inputs with dif- 
ferent values of ft it can be understood that good performance 
for the presented system can be obtained only for ft<ft 0 where 
ft 0 is the parameter belonging to that ellipse which has t lie same 
radius of curvature as the curve through (y,y') = (1,0). 
(Note that the radius of curvature of the approximating 
curves jumps at y' = 0 and is smaller for i/'< 0.) In the 
first quadrant (i. e., approaching y' = 0 from y'>0) this 
radius of curvature is given by 

P max == |p|j,' = +0=7*2/ 

Since at this point the ellipse radius of curvature is ft 2 /l, 

^ 2 <7nl7 

Therefore, for the first quadrant 

0<Vy3 (23) 

Here it should be noted that this is the same result as that 
obtained in the section entitled “Performance Evaluation of 
Sinusoidal Inputs” even though the parameter values are 
different. For the system of figure 20, y 3 =0.5; thus ft<CVb-5. 

Another example is given by the input x=l—e~ aT with 
d(x')/dx= — a. Since the smallest angular sector is at i/ = 0, 
a=2D(l-\-p max ) determines the limits on performance. 
This means that for the system represented in figure 20 good 
control can be expected for a value of a slightly smaller than 
0.75. 

A step input is represented by x= l—e~ aT with a— > 00 . The 
picture in the phase plane is x=l — (l/a)x' with a— > 00 . For 
very large values of a this is a straight line which forms a 
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Figure 20. —Superposed input and output phase planes showing avail- 
able angular sectors for sinusoidal inputs. i/3 = 1.5; ,7 = 0.45; 2 0 = 
0.5; 2T = 0.05; 2D = 0.25; co„=2; |0| maz = 2; \y\ max =0.5; |/3| min =l; 


small angle with x' = 0. For a-^co the curve degenerates to a 
point (x= 1, £' = 0). There is no doubt that a perfect followup 
of a step is not possible because for practical reasons the line 
for large values of a does not lie in the allowed angular sector 
at any point it is passing through. 

Related to the step input is the square-wave input. A por- 
tion of the output phase trajectory for a square-wave input 
to the system of figure 8 (e) is shown in figure 21. Figure 21 
is computed with the help of the differential equation 

<iy = y\ 

dy\ 2Z?(1 +W ?/ 1 + (1 +7 n )y-x 

This equation is obtained from equation (8) by replacing 
(Py/dr 2 by y n (dy\/dy ) . In this example one cannot immediately 
use approximating curves for designing the output because 
the error is too big at the start of the motion. The computed 
diagram is in good agreement with the test run shown in 
figure 8 (e). 



Figure 21. Portion of output phase diagram for square-wave input. 

0 3 = — /3 0 =2; 0 2 =— 0i = O.5; 73= — 7o=2; 72= —71 = t). 5; D=0.6; 

Cd, = 1. 

EFFECTS OF SWITCHING IMPERFECTIONS 

In equation (8) it is assumed that the parameters and y n 
change upon exact zero coincidence of any one or more of the 
variables ?/, y' , e, and e' as defined in equation (13) or by the 
logic of table 2. Physical imperfections, however, preclude 
this possibility. Thus, in the simulated model of equation 


(8), exact zero coincidence cannot be detected because of 
threshold effects, and relays are subject to time lag (both 
mechanical and electrical), dead zone, and chatter. 

In the present section an attempt is made to evaluate the 
effects that these switching imperfections had on the experi- 
mental result in order to obtain some practical design criteria 
for specifying switching requirements for good performance 
in other applications of this method of nonlinear control. 

The experimental results of figures 6, 9, and 11 can be used 
to demonstrate that relay imperfections were not important 
in the simulated model. Comparison of the sinusoidal re- 
sponses in these figures has shown that the period of the error 
became smaller as the magnitude of the input amplitude was 
increased. (This was discussed in the section “Performance 
Evaluation of Sinusoidal Inputs.”) From figure 11 (a) in the 
region of \y \ ^20 volts (the largest magnitude considered in all 
the experimental studies) the period of the error T e is 
measured as 

y # =«A«o.2 

There are four parameter switchings per error cycle (see fig. 
18). Assuming these to be approximately uniformly spaced, 
the minimum time between parameter switchings is approxi- 
mately TV4. Now assuming that the relays must be capable 
of closure in at least one-fifth of this time, the maximum 
allowable (real time) switching delay t md is 

u v tm& — T e /20=10~ 2 (23) 

In this simulated model 1 radian per second. Therefore, 
tmd — 10 millisec 

As given in table 4, the relays employed were capable of 
closure in 3 milliseconds or less so that they were entirely 
adequate for the experimental studies. 

(After the investigation reported here in the main text was 
finished, Dr. H. E. Lindberg performed a number of experi- 
ments to study quantitatively the influence of time delays on 
the performance of this control system (see appendix A). He 
investigated the influence of switching delays in the indi- 
vidual relays and additionally that of a delay in the output 
line of all relays (represented by an analog-computer approxi- 
mation of e~ t d p ). The results indicate that for sinusoidal 
inputs with frequencies up to 0.8 of the natural frequency of 
the undamped linear member co„, for example, the maximum 
error can be held within 5 percent of the maximum value of 
the input with relay delays as large as 47/co„ milliseconds. 
For a natural frequency of 10 cps this would be a delay of 
0.75 millisecond. After t his relay delay is selected, the maxi- 
mum error per input cycle is considerably less for lower input 
frequencies. For instance, at co=0.1co v the maximum error 
per input cycle would be slightly more than 1 percent. Also, 
the average error would be less than this (see fig. 31).) 

The ability to disregard the relays in the evaluation of the 
effects of switching imperfection on performance leaves only 
threshold effect in sensing the sign of the variables y, y\ e, 
and e f to be considered. As has been observed and discussed 
in the section “Performance Evaluation of Sinusoidal Inputs” 
the amplitude of the error for the nonlinear system was rela- 
tively independent of the magnitude of the input (i. e., when 
the system is operating in the rapid e to e' switching sequence 
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so that this is the minimum-error case). This constancy of 
the lower limit on the error magnitude is caused by switch- 
ing imperfections and thus threshold in sensing the signs of 
variables. Since under a normal switching sequence there 
are many more e and e' switchings than y or y' switchings 
and since the error is of an order of magnitude less than the 
error derivative, the primary cause of the lower limit on 
error is localized as threshold in sensing the sign of the error. 
Figure 13 shows the results of an investigation of this 
threshold effect. In all system-response figures except 
figure 13 the peak-to-peak threshold was approximately 14 
millivolts. In figure 13 (a) the peak-to-peak threshold was 
44 millivolts so that figures 7 (a) and 13 give a good com- 
parison of the effects of these two threshold values. 

CHOICE OF PARAMETER VALUES 

The performance of a completely specified nonlinear sys- 
tem has been discussed. The parameter values for this 
system, that is, 

D= 0.6 

03= — 00=2 

02= 0i = 0.0 

73 =— To=2 

72=— 71=0.5 

were initially chosen in the following manner: V was first 
selected to give good linear system performance; then the 
smallest values of the 0 m and y n parameters giving good non- 
linear system performance were chosen experimentally from 
a systematic variation of parameters utilizing the simulated 
system. This particular set of and y n parameter values 
can thus be denoted as ( 0 W , 7 n ) mfn since they establish the 
lower bound on parameter values for good nonlinear system 
performance. The physical significance of ( 0 m ,7 n )min is that 
loop gains and acceleration requirements of the linear mem- 
ber are minimized since 0 m and y n are feedback gain con- 
stants. 

From the discussion in the section “Performance Evalua- 
tion of Sinusoidal Inputs” centering about equation (18) or 
from the phase-plane methods of the preceding section it is 
to be expected that a general increase in the parameter 
values over ( 0 m ,7 n )min will result in improved nonlinear sys- 
tem performance by increasing the available acceleration of 
the system or increasing the angular sectors in the phase 
plane. Inspection of figures 16 (a), system 1, and 16 (c) y 
system 5, shows this to be true. In figure 16 (a) the param- 
eter values were ( 0 m ,7 n )min as given above, and in figure 
16 (c), system 5, they were 

03= -0o= 10 

02 = - 01=1 

73=— 7o=10 

72=— 7i=l 

On comparing the performance of these two systems it is 
noted that the corners of the input triangular wave are re- 
produced with less error by the system of figure 16 (c). 
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Thus far then it would seem that there is no upper bound 
on the parameter values; that is 

(0w>7n)maar ^ 00 

Invariably, however, there will be upper bounds on the values 
of the parameters because of acceleration limits in the phys- 
ical system. Figure 14 shows the effects of placing pro- 
gressively smaller acceleration limits on the system. Here 
it is seen that for |y"|z,>0.3|?/"| wax performance is not ap- 
preciably affected but for values less than this good per- 
formance is no longer obtained so that acceleration limits 
definitely tend to determine ( 0 W ,7 n )max- In general, then, 
there will be a whole range of values of p m and y n ; that is, 

(0W)7«)maar/ > (0wj7re) ^ (0 m>7 «) min 

for which good nonlinear system performance results. The 
final choice must depend upon the particular application and 
can easily be found experimentally. 

There are certain special cases of the 0 W and y n parameter 
values that lead to simplified feedback switching circuits 
and thus lead to nonlinear systems that are simpler to realize. 
In table 3 three of these are denoted as case (1) no derivative 
feedback, case (2) no e r sensing, and case (3) no e sensing 
in derivative feedback loop and no e' sensing in proportional 
feedback loop. 

By making inoperative the appropriate relays in figure 5 
the simplified switching circuits for these cases are easily 
obtained. It is desirable then to inspect the performance 
of these special cases (in comparison with that of the com- 
plete system) to ascertain whether or not as good performance 
can be obtained with less complexity. Figure 16 shows the 
response of these cases in comparison with that of the com- 
plete system. Figures 16 (a) and 16 (b) differ from figures 
16 (c) and 16 (d) in that the general magnitude of the 
parameters was increased in figures 16 (c) and 16 (d). 

These results show that case (2) is not worthy of much 
consideration since the magnitude of the oscillating error 
is large. Cases (1) and (3), however, should be considered 
for certain applications. For example, if it were known 
that the amplitude distribution of the input was relatively 
void near zero, then case (1) would serve as well as the 
case of the complete system. Case (3) shows nearly constant 
percentage error so that it could be useful in cases where 
accuracy was not so important as economy in components. 

To tliis point the parameter D has received little attention 
mainly because its value (within limits) is not particularly 
important. It lias been observed experimentally that D 
may be anywhere in the range O^W^l even for ( 0 w , 7 w ) m <n 
and performance of the nonlinear system is not affected. 
Inspection of the block diagram of the nonlinear s}^stem 
(fig. 4 or 5) shows that the physical significance of the 
variation of D for given values of 0 m , 7 * is that the damping 
factor of the linear member and the derivative feedback 
around the linear member change simultaneously. A case 
of more practical importance such as might arise in aero- 
dynamic applications of this type of control system is the 
variation of D for given values of 7)0 m , y n . That is, the 
damping factor of the linear member alone varies while the 
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feedback values around this member remain unaffected. 
Even in this case it lias been found that the nonlinear system 
performs well. For example, figure 17 shows the response 
of the nonlinear system 


D= 0.1 
ft — — ft=12 
ft=-ft=3 
7s=— 7o = 12 
72 =— 7 o = 0.5 

to triangular-wave inputs, while figure 7 shows the response 
of the system 


77 = 0.6 
ft = — ft =2 
ft = — ft =0.5 
73=— 7i = 2 
73=— 7 o =0.5 

to the same inputs. Comparing these results it is seen that 
performance is nearly the same for both systems. 

PRELIMINARY EXTENSION TO HIGHER ORDER SYSTEMS 

From a practical standpoint, limitations in the appli- 
cability of the nonlinear control system described by equation 
(8) do not stem from the inability to realize the feedback 
switching circuit but rather from assuming the linear 
member to be of second order. In many cases a more 
realistic approach is to consider the linear member to be of 
higher order but still predominantly second order. Figure 
22 shows the block diagram of a third-order system of this 
nature. The linear feedforward member could be, for 
example, the servo of figure 1 including the effects of arma- 
ture inductance. The same second-order feedback switching 
circuit was still employed. However, it could not be 



Figure 22. Block diagram of third-order nonlinear control system. 




(a) Nonlinear system. /? 3 = — / 3 0 = 2 ; / 3 2 = — Pi = 0A ; 73= — 7o=2; 72 = 

7i = 0.5. 

(b) Linear system. 

Figure 23.— Linear and nonlinear third-order system responses to ran- 
dom input. T=0.1; Z) = 0.6; 2.5 small divisions on time scale=l 
normalized time unit. 

expected that the P, n ,y n parameter values remain the same. 
Figure 23 compares the response of a third-order linear and 
nonlinear system. Here it is seen that the nonlinear system 
still responded with much less error than the linear system. 

On the basis of the comparative performance of these two 
systems it seems important to further studies toward control 
of general higher order systems using techniques similar to 
those developed in this investigation. 

CONCLUDING REMARKS 

From control equation (8) a second-order nonlinear 
control system that tends to maintain small instantaneous 
error for relatively arbitrary inputs has been synthesized 
using digital-computer techniques. The only type of 
nonlinear operations required in the realization are switching- 
type operations (zero-coincidence detection and parameter 
switching). The switching requirements are severe from 
the aspect of detector sensitivity and switching time delay 
but not impractical since all the switching is done in feedback 
paths at low power levels and may thus be performed 
electronically. 

The system demands sensing of error and error rate of 
change. Since general noise in a system of this type has to 
be expected, both error and error rate are smoothed. The 
influence of time delay in the relays and of a general transport 
delay on the performance of the system has been investigated 
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and diagrams have been developed which allow one to 
estimate the maximum error due to imperfections. 

It is a difficult task to compare the performance of a non- 
linear system with that of a linear system because no general 
criterion for comparison is available. Since for nonlinear 
systems the law of superposition does not hold, it is not 
adequate to choose the response to a certain input (e. g., the 
step input) as a criterion for performance comparisons. A 
number of different inputs therefore have been chosen for 
exhibiting the qualities of the nonlinear system. 


Experimental results indicate that this type of nonlinear 
control system performs better than a linear control system 
having a normalized frequency 15 times greater. Perform- 
ance is evaluated in terms of the average value of the magni- 
tude of the instantaneous error for band-limited inputs. 
Further, the nonlinear system performance is virtually 
independent of variations in the damping factor of the 
system. 

Stanford University, 

Stanford, Calif., October 18 , 1957. 


APPENDIX A 

INFLUENCE OF TIME DELAYS 

By H. E. Lindberg 


In tlie experiments reported in the main text of this report, 
the relays were carefully selected in order to keep the bad 
effects of finite closure time to a minimum and no really^ 
quantitative study of such closure delays was made at that 
time. Later, studies of this particular system were continued 
and quantitative information about the influence of imper- 
fections on the performance of the control system is presented 
here. 

One numerical and two experimental methods of finding 
response with relay delays are used. Plots of maximum 
error incurred per input cycle are given for sinusoidal inputs 
of various frequencies and for various relay delays. Effects 
of a transport time delay are also discussed. 

INFLUENCE ZONES AND TYPES OF DELAYS 

The system described by equation (8) can be interpreted as 
describing either of the s}~stems shown in figure 24 if there 
are no imperfections to be considered. However, there are 
certainly relay imperfections and there may be transport 
delays. Relay imperfections would affect directly those 



''Switching computer 



'S witching computer 

(a) Motor-load system. 

(b) Simplified airframe system. 

Figure 24. Block diagrams of two physical systems giving same 
differential equation providing no imperfections exist, m, point 
transport delay is introduced into circuit. 
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parts of the systems that are encircled by dotted lines in 
figure 24, that is, only the discontinuous feedback paths. 
If no further imperfections are considered, the results of such 
a study apply equally well to both configurations shown in 
figure 24 because the switching computer sees the same 
linear member in either case. Other studies, however, can 
conceivably be made on these systems where the difference 
between the two configurations becomes significant. A 
study of transport delays introduced at points m, for instance r 
would certainly give different results for figures 24 (a) and 
24 (b) because in figure 24 (a) the output y is fed directly 
through this point while in figure 24 (b) no such feedback 
is present because the equivalent of such a feedback is built 
into the differential equation of the mechanism being 
controlled (equivalent, i. e., when there are no imperfections). 

One of the most serious imperfections to be considered in 
the design of a relay servomechanism is the time lag between 
the application of a signal to the relay and the actual closing 
or opening of the relay. Associated very closely with this 
imperfection are other phenomena, such as signal threshold, 
relay dead time (the period of time in which the relay arm 
is in contact with neither terminal), and contact bounce. 
The only imperfection studied here is relay time lag, because 
in most relays dead time is very small as compared with the 
relay lag time. Contact bounce is also a secondary effect 
and is better studied in an actual installation. By appro- 
priate interpretation, the results of this time-lag study can 
be extended to include the effect of signal threshold which is, 
of course, a variable time delay. 

Rather than incre the natural frequency of the 

system being studied until the effect of inherent relay lag 
became noticeable, an artificial lag was introduced while 
the system was allowed to operate in a conveniently low 
frequency range. Two techniques were used to provide 
the time lag. The first was a direct delay of each relay 
signal, effected by providing an additional relay in series with 
each original relay and a resistance-capacitance delay circuit 
in between (this will be referred to as the electromechanical 
simulation). The second method consisted of inserting an 
analog-computer approximation of e~ tdV in the total output 
line of all of the relays. The two methods gave consistent 
results. 

ELECTROMECHANICAL SIMULATION OF TIME DELAY 

The armature of each logic relay in the control system 
with no intentional delays was energized by the plate current 
of a zero-coincidence detector circuit. In the delayed system, 
the plate current is made to flow through an added relay. 
This relay energizes an adjustable resistance-capacitance 
time-delay circuit which in turn feeds the armature of the 
logic relay. Two of the four circuits are shown in figure 25. 
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The amount of delay between activation of the zero-coin- 
cidence circuit and the closing or opening of the logic relay 
was measured by observing the coincidence pulse and the 
logic-relay contact pulse with a long-persistence oscilloscope 
whose beam traversed the scope at a known rate and was 
triggered by the coincidence pulse. The voltage supply in 
the time-delay circuit was adjusted so that the closing and 
opening delays were roughly equal. It is estimated that the 
accuracy of time-delay measurement was of the order of 10 
percent. Delays of less than 25 milliseconds ( t=r with the 
analog-computer wiring used) were difficult to measure by 
this method, and, more important, could not be easily simu- 
lated by this particular circuit because of high contact cur- 
rent in the added relay which caused the contacts to weld. 

The delayed system was operated using sinusoidal inputs 
of 20-volt peak-to-peak amplitude at frequencies from 1/10 
the natural frequency of the undamped linear system to the 
full value of the natural frequency. Nondimensional time 
delays of co+*=0.05, 0.075, 0.100, and 0.150 were used. Know- 
ing that the error caused by these delays should vary linearly 
with the applied voltage (see the section “Response to sinus- 
oidal inputs’ ’ which follows), runs were also made at 10- and 
40- volt peak-to-peak amplitudes as a check of the self- 
consistency of the data. 

Analog-computer simulation of time delay. — An analog- 
computer approximation of e~ l dV was wired using a circuit 


suggested by Morrill (ref. 6). The circuit is based on the 
Pade second-order approximation of e~ l dP 
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t*p 7 i 
1,728' 


. . (Ala) 


which is roughly equivalent to the first six terms of the Taylor 
expansion of e~ l d p about t d p=0: 




uy 

2 


t d f 

6 


uy uy 

24 120 


uy 

720 ' 


(Alb) 


The circuit is shown in figure 26. The element values were 
found by comparing the transfer functions of the circuit with 
the desired Pade polynomial fraction and equating coeffi- 
cients of powers of p. The transfer function of the circuit in 
a nondimensional form is: 


00 10 
0i 


Rio R4R10 R4R10C1 

Ri R3R7 R 7 p 


r 8 r 5 r 4 r 8 c 2 
R 2 Rc i ~ R 3 


P+(R 4 R s C 1 C 2 )p 2 


(A2) 


where the resistances R n arc as shown in figure 26. 
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Figure 26. — Analog-computer circuit for Pad6 second-order approxi- 
mation for e~ l d p . 


Comparing this with the Pade approximation written 


12 ^ 

12+6 tap+t/p 2 


(A3) 


6 equations for the 10 resistors and 2 capacitors can be writ- 
ten in terms of the time delay t d . The resulting six degrees 
of freedom on the circuit parameters were used to select the 
convenient values shown. Figure 27 shows the location of 
this delay circuit in the overall simulation of the relay con- 
trol system. 



Figure 27. -Insertion of time-delay circuit of figure 26 into control- 

system analog. 

Since this method was used only as a check on the electro- 
mechanical method, the only input studied was a constant 
of 10 volts at delays of t d = 0.025, 0.050, 0.075, and 0.100 second. 

Simulation of a time delay by this method allowed greater 
flexibility of the location of the time delay. The delay was 
inserted as shown in figure 28 to simulate the effect of a 
transport time delay in the system. The output for this 
configuration was essentially the same as that found for relay 
delays. This result was to be expected because the essential 
difference between a transport delay and relay delay in this 
system is that relay delays affect only the feedback terms 
which go through the relays while transport delays affect 
the input x as well. But for the great majority of operation, 
x changes very slowly as compared with the changes in the 
feedback terms and lienee, it makes very little difference in 
the response whether the input x is delayed. 



Figure 28. — -Insertion of time-delay circuit for study of transport 

delays. 

The responses of the two configurations become more 
nearly identical as the magnitude of the time delay decreases. 
If the basic controlled element is a motor load as shown in 
figure 24 (a) and if a transport delay were introduced at 
point m, the resulting output would differ somewhat more 
from a relay-delayed output because then both x and the 
direct y path as well as the discontinuous relay signals would 
be delayed. No tests were run with this configuration. 

Computational study of relay delay. -Response of the 
system with relay lag was computed using the differential 
equations directly. To simplify the amount of calculations, 
only response to constant inputs were computed. Calcula- 
tions were based on approximating the output with a series 
of parabolic arcs. This is equivalent to assuming that y n 
is constant during each computation interval, a reasonable 
assumption except for large time delays. The input studied 
was #=10 and calculations were made until the limit cycles 
dE 

on a plot of E against were approached within the accuracy 

of computation. Computations were made such that the 
limit cycles were approached from the inside in some cases 
and from the outside in others. Limit cycles for various 
time delays are shown in figure 29. Computation intervals 
were such that two parabolic arcs were used in each quadrant. 

The computation details are as follows: With the assump- 
tion that the output is a series of parabolic arcs, the equation 
for the error during the nth interval between switch points is: 

En=E ni +E ni 'T n +tE ni ''T* (A4) 

where r n = 0 at the start of any interval and E nt , E n /, and 
E ni " are the initial error, error velocity, and error accelera- 
tion, respectively, for the nth interval. Upon differentiating, 
the velocity during any interval becomes: 

E n '=E ni '+E ni "r n (A5) 

Eliminating the time r n between these two equations, the 
equation of the phase-plane trajectories becomes: 

tt r rr __I (jr. ')2 ( Ac) 

L ni 2 EJ J 2 E ni " [ n } 1 } 

Equations (A4) and (A5) are used to find E as and EJ 
at an actual switch point by letting T n =u v t d (the delay time) 
and using as initial conditions E at E' = 0 or E r at E=0 
(the conditions at a theoretical switch point). When this 
point is found, equation (A6) is used with E ni =E as and 
E n { = EJ and the next theoretical switch point is computed 
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Figure 29. — Error limit cycles due to switching delays for a constant 
input of £=10 volts. 


by setting either E n = 0 or E n '= 0 depending on which type 
of switch point is being approached. In all cases the accel- 
eration E ni " is found from the differential equation of the 
control system, equation (8). For the special case where 
£=10, equation (8) can be rewritten using y=E-\-x as: 

E"+2D(l+p m )E' + (l+y n )E=-10y n (A7) 

An effort was made to determine some general method for 
predicting the limit cycles as an analytic function of t d , D , 
j3 m , and y n but, even with this very simple input, the 
limit cycles are approached in such a complicated way that 
finding such a function would be a formidable task at best. 
However, a systematic series of calculations as shown here 
could be made and charts plotted to indicate the effect of 
these various parameters. 

Comparison of the three methods for constant inputs. 

For a constant input of £=10 a good comparison of results 
from the three methods of studying relay delay was made by 
constructing a logarithmic plot of error versus time delay. 
Data from the two experimental methods had to be modified 
slightly to account for actual relays delays and thresholds 
of zero coincidence. The modification consisted of measur- 
ing the system error with no intentional delay, finding the 
delay required to produce this error from the calculated plot 
of error versus delay, and adding this delay to each of the 
experimental points. (The effective unintentional relay 
delay was found to be 20 milliseconds. This is considerably 


greater than the 3-millisecond delay specified by the relay 
manufacturer. The additional effective delay is due to the 
inherent thresholds in the switching circuit which can be 
considered as a variable time delay. Also, the actual delay 
of the relays may have exceeded 3 milliseconds.) 

After this modification, the analog-computer data were 
indistinguishable from the calculated data. Data from the 
electromechanical method also agreed quite well but differed 
from the computed results by as much as 10 percent at a 
time delay of 0.12. (See fig. 30.) 

Response to sinusoidal inputs. As mentioned before, the 
error resulting from relay lags varies linearly with the magni- 
tude of the input x. That this is true is easily seen by examin- 
ing equation (8). If x is replaced by Kx and y by Ky where K 
is a constant, the differential equation remains unchanged 
and the error e=x—y is now e*=Kx—Ky=Ke. 

The reader should be warned that although this is a well- 
known property of linear differential equations, this system 
remains essentially nonlinear. Using the definitions of 
and y n given in the text following it, equation (8) can be 
written : 

y"+2Z)(sgn y'— & sgn e— 2 & sgn e')\y'\ + 

(sgn y — jy sgn e— 2 y sgn e')\y\=x (A8) 

With the equation in this form, the reader can easily verify 
that if y x and y 2 are solutions to equation (A8) in its homoge- 
neous form (£=0) then y x -\-y 2 is not necessarily a solution 
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Figure 30. — Comparison of methods of switching-delay study with a 
constant input of £=10 volts. 
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(a) co = 0. lo>^. 

(b) Cl) = OAUy, 

(c) co — 0. Goo,,. 

(d) co = con- 
figure 31.— Responses for inputs of various frequencies with relay delays of co^d = 0.10. 


because of the nonlinear character of the absolute value and 
signum functions. 

However, it can be concluded from the linear property 
that if the error encountered for any constant input is 
obtained, the error for any other constant input can be found 
by a simple scale change. For very low frequency sinusoidal 
inputs where x' and x " are relatively small it was found that 
the error was governed almost entirely by this consideration. 
This is demonstrated by figure 31 (a) for a 40-volt peak-to- 
peak input with co = 0.1a>„ and a delay of co„£ d =0.100. Notice 
that the envelope of the error is very nearly a sine wave with 
its peak value occurring at the peak value of the input. 

As the input frequency is increased, the point of greatest 
error shifts away from the peak of the input. In figure 31 (b) 
for co=0.4a>„ the greatest error occurs almost % cycle after 
the input reaches its peak in either direction. This is because 
the (1 term becomes sizable in this frequency range 

during the acceleration period (i. e., as \x'\ increases) and 
reaches a maximum at l / 4 cycle past the input peak. This 


large feedback term is applied in an undesired direction for 
0.10 unit of nondimensional time and causes the large lag 
error. 

At still higher frequencies the y" term coupled with the 
(l+ 7 „)y term both reach their maximum at the peak value 
of the input and the greatest error again occurs near this 
time. An example of this is given in figure 31 (d) for an 
input frequency of w=1.0oo„. For intermediate frequencies 
the peak shifts around depending on initial conditions as 
shown for co=0.6a>,, in figure 31 (c). 

From an overall standpoint a plot of the greatest error 
versus frequency showed an approximately linear increase 
of error with frequency (see fig. 32). Similarly, a plot of 
greatest error versus time delay is linear up to about 
co„/ d =0.10 for frequencies up to co=0.8co„ (see fig. 33). An 
interesting comparison of these results and the results found 
for the system with no intentional delays is shown in figure 
34. Also plotted in this figure is the maximum error for the 
linear system alone (y" -\-2Dy' -\-y=x). 
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0 .1 2 .3 .4 .5 .6 .7 8 .9 1.0 


Frequency, oj/w v 

Figure 32. — Maximum error occurring during each input cycle plotted 
against input frequency with switching delay as a parameter. 

rc=10 sin cot = 10 sin^— volts; /S 3 = — /3 0 = 2 ; 

02= — 0i = O.5; 73= —72=2; 

72= —7i = 0.5 ; Z) = 0.6; co „= 1.0. 


Specification of time-delay tolerances. — If an actual sys- 
tem were to be set up such as that described by equation (8), 
figure 32 would provide an indication to the designer of how 
to specify relay-delay tolerances. Take, for example, a 
system which has an undamped natural frequency f v = 10 cps. 
If the error of this system is to remain below 5 percent of the 
maximum value of the input for frequencies up to co=0.8co„, 
then the relay delays must be held below 0.75 millisecond. 
This result is found by first finding that value of co v t d on the 
co=0.8co„ curve of figure 33 at which the maximum error is 
5 percent. This value of c o v t d is 0.047. The actual time delay 
t d is: 


0.04/ 0.04/ 0.04/ n nnn7K 

t d = = ~o~~r~^on — =0.00075 sec (A9) 

to* Zi rfy Z U7r 


This value of time delay will not necessarily assure that 
the error will always remain within 5 percent of the input 


for other types of inputs which are band-limited to 0.8a>, 
because of the nonlinear nature of the system. However, 
it is known from the special type of inputs studied that if 
the error is to be as specified the time delay can be at most 
0.75 millisecond. 

For sinusoidal inputs with co<0.8co„ the maximum error 
will be less than 5 percent of the maximum value of the 
input. For example, at w^O.lco, the maximum error will 
be slightly more than 1 percent . 



0 0.05 0.10 0.15 

Relay delay, w v t d 


Figure 33. — Maximum error occurring during each input cycle plotted 
against switching delay with input frequency as a parameter (same 
data as for fig. 32). 



Figure 34. — Maximum error plotted against input frequency for 
several time delays and for linear member only. 


APPENDIX B 


COMPARISON OF EQUATION (8) AND NONLINEAR CONTROL EQUATION OF REFERENCES 1 AND 3 TO 5 

COMPARISON OF NOTATIONS 


In original notation the nonlinear control equation was 
given as (refs. 1 and 3 to 5) 


a+ ^+ b ji+ c yo=y<w 


(Bl) 


&=& + £l-fp 1 sgn (y 0 'E)+^ sgn (*//£') J 
c=c + 1~ 


1+^F sgn (y 0 E)+~ sgn 


E=(y 0 —Vi) 

( )'=<*( )/dt 


where 

Vo(t) 

y & ) 


output 

input 


Abi,Ab 2 ,Aci,Aco constants, Ab 2 ,Aci^> Ac* 

a + ,6 + ,c + constants 


30 


REPORT 1391 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


In terms of the notation of equation (8) the above equation 
is written 


® + ^+& + (l+&») ^+c + (l+y n )y=x(t) (B2) 


where 


0m= — I0sgn (y'e)— 2 p sgn (y'e');m= 0, 1, 2, 3 
y n = — iTSgn (ye)— 2 y sgn (^');/i=0, 1,2,3 
e =(x—y) 

i0, 2/3,17,27 constants, ij8> 2 0 , i 7>27 


a + ,b + ,c + constants 

The subscript convention is 


03 = 10 + 20 

73 = 17 + 27 

02 = 10 20 

72 = 17 — 27 

01 = -02 

7i=— 72 

00= 03 

e 

II 

1 

Evidently, the correlation betw 
is then 

een equations 


(Bl) and (B2) 


V=y< 


In either notation a set of seven (constant) parameters is 
needed to characterize the system. For example, in equation 
(B2) the set (a + ,b + ,c + , i0, 2 0,iy, 2 y) is sufficient. 

NORMALIZATION OF NONLINEAR CONTROL EQUATION 


When considering control systems it is possible to reduce 
the number of parameters necessary to specify the nonlinear 
system. To do this a form of normalization familiar to 
linear theory is utilized. First, it should be realized that if 
error is to be defined as e=(x—y) then c + = 1. Thus, 
assuming c + = 1 equation (B2) may be written 


where 


1 d 2 y . 2D 1 _ >. dy . ✓ . 

dt +il+yn)y=x{t) 




2 



CB3) 



0) v 


Introducing normalized time r=co„£, equation (B3) be- 
comes 

-j— -\-2D(l+p m ) j^+(l+y n )y=x(T) (B4) 


i0= 


A 6 , 


X = l/i 

a A b* A Ci A Co 

2/3_ b+ ’ l ' Y ~7 r ’ 2y ~'c+ 


Thus, knowing the natural frequency of the undamped 
linear system, a>„, the number of parameters necessary to 
specify performance in the nonlinear case is five, that is, 


e=-E 

and a + , b + , and c + are the same. 


or alternatively 


(D, 2 0, 17, 27) 

(D, 03, 02, 7s? 72 ) 


APPENDIX C 

SIMULATION OF SECOND-ORDER NONLINEAR CONTROL SYSTEM 


EQUIPMENT 

Experimental studies were carried out with the following 
equipment: 

(1) General equipment: 

(a) Analog computer, Beckman Ease 

(b) Low-frequency function generator, lip 202 A 

(input device) 

(c) Pen recorder, Sanborn Twin-Viso (output 

device) 

(d) Vacuum-tube voltmeter, RCA WV-97A 

(e) Oscilloscope, Dumont 304-A 

(2) Special equipment: 

(a) Binary-logic switching circuit employed in 
conjunction with (a), described in detail in 
section “Binary-Logic Relay Switching 
Circuit. ” 

COMPUTER SETUP 

The computer diagram for the differential equation 

^+2Z?(l+0 m ) ^+(1+7 n )y= x ( T ) (Cl) 

is given in figure 35. The correspondence between this sim- 
ulated system and the block diagram of figure 4 or 5 is 


straightforward (see ref. 7). Operational amplifiers 1 to 6 
are used in simulating the linear member of the physical 
nonlinear system. Resistor Rj provides adjustment of the 
linear damping factor D. The input to this simulated linear 
member is x— (2Dp m y f -\-y n y) where x is obtained from the 
input device, while the values of feedback 2 D@ m y f and 7 n y 
are obtained with resistors R 2 through R 9 (see table 5) 
connected to y' and y through a binary-logic relay switch- 
ing circuit derived in the following section. Amplifiers 7 
to 10 are sign changers. The four variables y, y' , e, and e' 
whose signs are to be sensed are made available as shown in 
the lower right of figure 35. 

BINARY-LOGIC RELAY SWITCHING CIRCUIT 

To complete the simulation of equation (8) it is neces- 
sary to implement the binary logic of table 2. Figure 5 
shows a relay switching circuit (designed on an “and/or” 
basis) that realizes the necessary logic. 

In order to preclude the possibility of time delay in switch- 
ing, “fast” relays have been employed. The average prop- 
erties of the 14 double-pole single-throw relays comprising 
the switching circuit are given in table 4. 

The necessary synchronism in relays connected by dashed 
lines (fig. 5) is obtained by series connection of field coils. 


INVESTIGATION OF A NONLINEAR CONTROL SYSTEM 


31 



Figure 35. — Computer diagram for simulation of equation (8). 


Into the four channels of field coils thus formed is read the 
four-digit binary logic of table 2 where now 

0=Function>0=No coil current 

1 =Function<0=Coil current 

This process of reading in the binary logic may be done by 
monitoring the variable (y, y' , e , e') with four separate zero- 
coincidence detectors or amplitude selectors (denoted CD 
in fig. 5 and described in detail in the following section). 
The output of these coincidence detectors then drives the 
respective channels of relay coils. 

SENSITIVE TRANSISTORIZED ZERO-COINCIDENCE DETECTOR 

One of the four identical amplitude selectors designed for 
use with the switching circuit of figure 5 is shown in figure 36. 
This circuit basically consists of a grounded emitter PNP 
transistor stage T x driving a pentode pulse amplifier T 2 . 
(It should be possible to utilize a transistor in place of the 
pentode if desired.) Four relay coils connected in series 
form the plate load of T 2 . Positive feedback /? 3 has been 
incorporated for regenerative switching. 


-7 5 v B+ (210 v regulated) 



Figure 36. — Transistorized zero-coincidence detector. 


The operation of this circuit is as follows: When the input 
voltage V becomes more negative than —50 millivolts, 
current starts to flow in the base-emitter region of T x 
This initiates collector-emitter current which effectively 
grounds the collector and thus the grid of T 2 . Plate cur- 
rent flows in T 2 so that the relays which form this plate load 
are switched to the up position (fig. 5). Positive feedback 
i? 3 has been incorporated to make the switching regenera- 
tive. Stages T x and T 2 will continue to conduct until V 
goes positive by 50 millivolts, at which time current is cut 
off in T x and thus T 2 ; the relays return to the normally 
closed positions (down in fig. 5). 

Refinements that have been included in this circuit are: 

(1) Adjustable positive feedback R 3 . This feedback gives 
control over the zero-sensing t hreshold of the detector. The 
greater the positive feedback the greater the threshold. 
This adjustment is desirable for studying the effects of 
threshold in sensing discussed in the section “Effects of 
Switching Imperfections.” The measured peak-to-peak 
threshold values obtainable with this arrangement were a 
maximum of 400 millivolts and a minimum of 60 millivolts. 
(In the computer setup it was necessary to amplify error e 
that formed the input to one of these detectors five times in 
order to bring the switching threshold down to 15 millivolts, 
peak to peak.) 

(2) Adjustable bias for emitter of Ti. This is necessary to 
compensate for the slight positive bias (approximately 
equal to 100 millivolts) given to the base of Tj by the positive 
feedback R 3 . 

(3) A clamping or clipping diode T 3 to protect the transis- 
tor r i\ from excessive base-emitter inverse voltages. Thus 
it can be seen that the input impedance of this coincidence 
detector is 15 kilohms since the base of T x is effectively always 
grounded. 


32 


REPORT 1391 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX D 


BEHAVIOR OF OUTPUT IF INPUT AND INPUT DERIVATIVE ARE SIMULTANEOUSLY SMALL 


Assume that the input is a curve which may be approxi- 
mated by 


for values r<^5. Since x'=2ar, both input and input 
derivative vanish with r— >0. The input phase curve is 

x =l (*')* 

The output depends on the initial conditions: 

(1) For y( 0) = 0 and y'( 0) = 0 


y=a \_h 


4 i 


r 6 + 


•••] 


.12 1 30 

(2) For 1/(0) =0 and y'(0) = e l 

2/=€i[r— Z)(l+/? m )T 2 4- . . .] 


(3) For y(0) — ( 2 and y'(0)=0 

y=e2 [ 1+ (^) T2+ - * •] 

In all cases reproduction of the input is not perfect very 
close to r=0; hourever, in cases (2) and (3) switch points 
may occur for rather small values of r. In case (1) no 
switch point close to t = 0 can be expected. This case will 
rarely occur; in most cases neither y( 0) nor y'( 0) will be 
zero. Then the output is a superposition of cases (2) and 
(3). In this event the error e=x—y is given by 

e= — e 2 —e i T+^oi-\-eiD(l-\-p m ) —6 2 • • -J 7 " 2 

It is evident that for small values of €i and e 2 the error grows 
with a. This can easily be seen in the results of figure 11. 
The input was x=A(l— sin fir); that means, near x=0> the 
input may be approximated by a parabola with a=l/2(fi 2 )A 
or the errors near a?=0 grow with fi 2 . 


APPENDIX E 

APPROXIMATING CURVES 


The approximating curves which can be used to trace the 
output for a given input form a network in the phase plane. 
It has proved practical to present the network for e and e' 
with the same sign in one sheet (e. g., see fig. 20) and the 
network for e and e ' with opposite signs on another sheet. 
Superposition of both sheets allows establishing of the 
four approximating curves through each point of the phase 
plane. 

The approximating phase curves are determined by 
equation (21). Since the values \ and X 2 change from 
quadrant to quadrant the curves are composed of portions 
of different analytical curves which are patched at y' = 0 
and y= 0. 

The roots and X 2 depend on D, /3 m , and y n ; they may 
be complex or real. If the roots are complex, the approxi- 
mating curves are curves of the spiral type which wind 
around the origin of the phase plane. In the case of real 
roots the curves have a quite different character. This 
can easily be understood by transforming the equation of 
the approximating curves (see refs. 1 and 3 to 5 and ch. V 
of ref. 8). For real roots the new coordinates n and v are 
introduced: y\-y'=y 

y*2—y ,= v 


Then equation (21) yields 

i^ 2 ==A/ju Xi 

If Xi and X 2 are real and of opposite sign, the curves in the 
\iv plane have a hyperbolic character with (n,v)= (0,0) as 
saddle point (see fig. 37 (a)). If and X 2 are of equal 
sign, the curves have a nodal point (see figs. 37 (b) and 
37 (c)). In figures 38 and 39 phase curves with saddle 
point and stable and unstable node are shown in the original 
y y f plane. 

The approximating curves are composed of portions of 
these different types (see fig. 40). 

There is no need of avoiding p m ,y n combinations which 
lead to node- type approximating curves because onty por- 
tions of these curves are used. In the earlier publications it 
appeared as if (in either one of the networks of approximating 
curves) one set of curves had to be formed by curves of 
spiral character (complex roots X). However, this has 
proved to be an unnecessary restriction. There might be 
some trouble with node-type curves if large delays in switch- 
ing should occur (e. g., delays in y' = 0 switching in fig. 40). 
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(a) 


(b) 





(a) £2 <0. 

Ai 

(b) r! > L 

(c)0< ^ <1. 

Figure 37. — Curves in nv plane. 





(a) Stable node. 

(b) Unstable node. 

Figure 39. — Phase curves with nodal point. 



Figure 40. — Composition of an approximating curve. 
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TABLE l — CODING SCHEME FOR SWITCHING LOGIC FROM 
EQUATION (14) 


[0, function>0; 1, function<0] 


Binary coded decimal 

Decimal 

y 

y ' 

e 

e' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

2 

0 

0 

1 

1 

3 

0 

1 

0 

0 

4 

0 

1 

0 

1 

5 

0 

1 

1 

0 

6 

0 

1 

1 

1 

7 

i 

0 

0 

0 

8 

i 

0 

0 

1 

9 

i 

0 

1 

0 

10 

i 

0 

1 

1 

1 1 

i 

1 

0 

0 

12 

i 

1 

0 

1 

13 

i 

1 

1 

0 

14 

i 

1 

1 

1 

15 


TABLE 2.— MATRIX OF ALLOWED 0 m , 7n COMBINATIONS 
AS DETERMINED BY ENCODED SWITCHING LOGIC 
OBTAINED FROM EQUATION (14)“ 



00 

01 

02 

03 

7o 

0000 

1111 



0100 

1011 

7l 


0001 

1110 

0101 

1010 


72 


0110 

1001 

0010 

1101 


73 

0111 

1000 



001 1 
1100 


a Examples: 

(1) /3370 is chosen when 0100 occurs, i. e., (i/>0, ?/'<(), e>0, e'>0), 

or when 1011 occurs, i. e., (t/< 0, */'> 0, e<0, c'<0); 

(2) 0 2 7i is chosen when 0101 occurs, i. e., (?/>0, ?/'< 0, e>0, e'<0), 

or when 1010 occurs, i. e., (z/<0, */'>0, e<0, e'>0); 

(3) 0i 73 is not possible. 


Figs. 16(a) and 16(b) 

Figs. 16(c) and 16(d) 

Complete system 

System 1 : 

03 — — 0o = 2 73 = — 7o = 2 

0i>= —01 = 0.5 72 =— 7i = 0.5 

System 5: 

03 =—0o=lO 73= — 70 = 10 

02 = — 01 = 1 72 = — 71 = 1 

Case ( 1 ) No derivative feedback 

System 2: 

03= —0o=O 73 = —7o=2 

02 = —0i = O 72 =— 71 = 0.5 

System 6 : 

03= — 0o=O 73 =— 7o=10 

02= —01 = 0 72= — 71= 1 

Case (2) No e' sensing 

System 3: 

03= —00 = 2 73= —70 = 2 

02= —01 = 2 72= —71 = 2 

System 7 : 

03= — 0o=lO 73 =— 7o=10 

02 = — 0i = lO 72 = — n = 10 

Case (3) Noe sensing in derivative feedback loop and no e f sens- 
ing in proportional feedback loop 

System 4: 

03= —00 = 2 73= —70 = 2 

02= —01= —2 72= —71=2 

System 8: 

03= — 0o=lO 73 =— 70 = 10 

02 = — 0 i= — 10 72 =— 7!= 10 


TABLE 4.— AVERAGE PROPERTIES OF RELAYS FOR 
SWITCHING CIRCUIT OF FIGURE 5 

Description: General Electric CR-2791 double-pole single-throw 

relay mounted on 5-prong Amphenol base with permanent alu- 
minum dust cover 


Connections: 




Electrical Properties 

Values 

Coil resistance 

2.5 kilohms 

Coil current for positive 

5 to 10 ina 

action 


Pull-in time 

2 millisec 

Drop-out time 

3 millisec 


TABLE 5.— CORRESPONDENCE BETWEEN ADJUSTABLE 
COMPUTER ELEMENTS AND PARAMETERS OF EQUA- 
TION (8) 


Element 

Corresponds 

to 

Units, 

kilohms 

Ri 

(2D)-' 

100 

R 2 

(2 DM-' 

100 

r 3 

( 2DM~' 

100 

1L 

(2 DM-' 

100 

Ro 

(2DM-' 

100 

Re 

( 70 ) 1 

100 

r 7 

(yi) 1 

100 

Its 

(72)- 1 

100 

R« 

(73)- 1 

100 
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